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Abstract 

Current  operational  needs  require  the  deployment  of  radiation  detection 
equipment  with  the  ability  to  accurately  and  reliably  identify  special  nuclear  materials 
and  their  byproducts  without  dependence  on  cryogenics.  This  requires  a  resolution  of 
0.5%  or  less  over  a  range  of  200  to  700  keV.  The  feasibility  of  a  Compton  spectrometer 
to  achieve  this  resolution  is  examined. 

The  Compton  spectrometer  system  used  consists  of  two  detectors.  The  Compton 
scatter  event  occurs  in  a  CdTe  detector  where  the  Compton  electron  energy  is  collected. 
Gamma  rays  scattered  out  of  the  CdTe  at  an  angle  detennined  by  a  conical  collimator,  are 
collected  in  a  Nal(Tl)  detector.  Coincidence  electronics  determine  correlated  events  and 
allow  the  Compton  electron  and  scattered  gamma  ray  energy  spectra  to  be  collected. 

Experimental  and  modeling  techniques  are  used  to  evaluate  the  system’s 
resolution  and  efficiency  and  provided  reasonable  agreement.  Expected  experimental 
results  based  on  previous  work  were  not  reproduced  and  the  source  of  the  difference 
remains  unknown.  Results  suggested  strict  requirements  of  collimation  will  make  some 
low  areal  count  rate  applications  impossible. 
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FEASIBILITY  ANALYSIS  OF  A  COMPTON  SPECTROMETER  SYSTEM  FOR 
IDENTIFICATION  OF  SPECIAL  NUCLEAR  MATERIAL 

I.  Introduction 

The  Department  of  Defense  (DoD)  has  many  activities,  such  as  nuclear  treaty  and 
counter-proliferation  monitoring,  that  require  deployable,  high-resolution  radiation 
detection  systems.  The  primary  focus  of  this  research  is  the  development  of  a  system 
capable  of  achieving  high  resolution  over  a  wide  range  of  energies  that  is  not  dependent 
on  resource  intensive  liquid  nitrogen.  The  system  analyzed  in  this  research  is  based  on 
the  previous  work  of  Captain  Williams  on  the  development  of  a  Compton  spectrometer 
system.  The  benchmark  used  to  evaluate  the  system’s  perfonnance  is  the  difficult  task  of 
determining  the  isotopic  ratio  of  a  plutonium  source,  as  needed  for  treaty  monitoring. 

This  application  can  easily  be  extended  to  other  DoD  radiation  detection  needs  requiring 
similar  energy  resolutions  over  a  wide  range  of  energies. 

Background 

As  the  face  of  world  politics  continues  to  change,  the  requirements  for  improved 
nuclear  counter-proliferation  and  treaty  verification  techniques  increases.  “Today,  the 
most  important  threat  to  US  security  is  the  proliferation  of  weapons  of  mass  destruction 
(WMD)  to. .  .countries  like  North  Korea,  Iran  and  Iraq  as  well  as  to  non-state  actions  such 
as  terrorist  and  organized  criminals”  (Hays,  1998:  2).  The  accountability  of  existing 
special  nuclear  material  (SNM)  and  the  detection  of  nuclear  technology  development  are 
critical  to  providing  security  to  our  nation  and  stability  to  the  world.  A  portable  detection 
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system  capable  of  identifying  the  presence  and  the  isotopic  makeup  of  SNM  is  needed  for 
both  of  these  missions. 

“On  May  24,  2002,  the  United  States  and  Russia  signed  the  Treaty  on  Strategic 
Offensive  Reductions  (also  known  as  the  Moscow  Treaty),  establishing  a  limit  of  1,700 
to  2,200  warheads  in  each  party’s  deployed  strategic  nuclear  arsenal  by  2012”  (Graham, 
2003:  1453).  The  accountability  of  weapon  systems  in  this  process  is  critical  as  they  are 
moved  to  storage  or  decommissioning  facilities.  With  this  overwhelming  burden  on 
accountability,  a  considerable  amount  of  time  may  elapse  before  a  missing  weapon 
system  is  noticed.  With  the  possibility  of  “spoof’  weapons,  replacing  a  weapon’s  physics 
package  with  non-weapons  grade  material,  the  missing  weapon  may  never  be  identified. 
The  ability  to  determine  the  isotopic  ratios  of  the  material  with  a  portable  system  would 
facilitate  Strategic  Arms  Reduction  Treaties  (START  I  &  II)  inspections  and  treaty 
verification. 

The  clandestine  nuclear  weapon  programs  of  countries  like  North  Korea  and  Iran 
pose  a  serious  threat  to  the  world  and  a  difficult  challenge  to  the  intelligence  community. 
The  development  of  un-accounted  nuclear  weapons  and  SNM  cannot  be  tolerated.  The 
detection  of  such  activities  can  be  determined  from  their  radioactive  byproducts.  Again, 
deployable  and  dependable  detection  systems  with  high  resolution  are  required. 

Problem  Statement 

The  grade  of  a  plutonium  sample  is  difficult  to  determine  because  of  the 
complicated  gamma  ray  energy  spectra  of  239Pu  to  240Pu  which  is  further  complicated  by 
the  effects  of  shielding.  The  grade  of  plutonium  is  classified  by  the  isotopic  ratio  of  Pu 
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to  240Pu.  Weapons  grade  plutonium  is  identified  as  having  over  90%  239Pu  and  less  that 
7%  240Pu  (DTRA,  2001:  120,  135).  The  high  spontaneous  fission  rate  of  240Pu  is 
undesirable  for  use  in  weapons  because  it  leads  to  early  initiation  and  reduced  yield  or 
device  failure.  Reactor  grade  plutonium  is  identified  as  having  15-25%  240Pu.  The 
complete  energy  spectra  of  240Pu  and  239Pu  are  found  in  Appendix  A.  Although  there  are 
many  characteristic  gamma  rays  of  both  isotopes,  the  639.99  keV  and  645.9  keV  ~  Pu 
peaks  and  the  642.35  keV  240Pu  peak  are  used  for  illustration.  The  spectra  of  239Pu  and 
240Pu  around  these  energies  are  shown  in  Figure  1  and  Figure  2,  respectively.  These 
energies  are  chosen  for  two  primary  reasons.  The  first  is  the  higher  probability  of  escape 
from  the  nuclear  weapon  casing.  Second  is  the  availability  of  peaks  from  both  isotopes 
in  a  narrow  energy  range.  This  limits  the  energy-dependent  effects  of  scattering  and 
absorption  (Williams,  2003:  5-6).  The  ability  to  distinguish  these  peaks  requires  a  full- 
width  at  half-maximum  (FWHM)  of  approximately  2  keV.  Although  this  resolution  can 
be  achieved  with  high  purity  germanium  detectors  (HPGe),  their  dependence  on  liquid 
nitrogen  (LN)  for  cooling  makes  them  less  practical  for  field  use. 
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Figure  1.  Energy  spectrum  for  Pu  showing  639  and  645  keV  peaks  of  interest 
and  the  complicated  spectrum  requiring  high  resolution  (INEE,  1999) 


Figure  2.  Energy  spectrum  for  240Pu  showing  642  keV  peak  of  interest  (INEE,  1999) 
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Scope 


The  goal  of  this  research  is  to  use  the  Compton  spectrometer  for  plutonium 
isotopic  ratio  determination,  continuing  the  work  of  Capt.  Chris  Williams  in  2003.  Capt. 
Williams  determined  that  the  Compton  spectrometer  showed  the  required  resolution  to 
separate  the  plutonium  peaks  of  interest,  see  Appendix  B.  His  research  used  an 

i  ■jn 

essentially  mono-energetic  Cs  source. 

The  purpose  of  this  research  is  to  investigate  the  resolution  and  efficiency  of  the 
Compton  spectrometer  system  over  an  energy  range  of  200-700  keV,  with  multi- 
energetic  sources.  The  ultimate  goal  is  to  reconstruct  an  operational  Compton 
spectrometer  system,  using  the  same  detectors  and  collimators,  with  two  major 
considerations  in  mind;  the  reproducibility  of  experimental  geometry  and  the  additional 
signal  processing  needed  for  multi-energetic  sources.  In  addition  to  the  experimental 
procedures,  a  simulation  code  was  developed  to  evaluate  experimental  results  and  to 
assist  with  design  modifications  for  follow  on  research. 

Investigative  Questions 

1 .  Can  the  previous  results  by  Capt.  Williams  be  reproduced? 

2.  How  does  the  energy  gating  affect  the  Compton  Spectrometer? 

3.  Can  modeling  the  Compton  spectrometer  help  with  future  design? 

Order  of  Presentation 

The  following  chapters  cover  my  research  of  the  Compton  spectrometer.  Chapter 
II  discusses  the  fundamental  physics  that  makes  the  system  possible  and  reviews  some 
essential  radiation  detection  measurement  theory.  The  equipment  used  during  this 
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research  is  described  in  Chapter  III.  Chapter  IV  covers  the  procedures  used  during  this 
research.  Chapter  V  presents  the  experimental  and  simulated  results  and  analyses.  The 
conclusions  from  the  data  analyses  and  recommendations  for  future  research  are  covered 
in  Chapter  VI.  Finally,  Appendix  A  through  P  present  detailed  background  infonnation 
including  tabulated  data,  computations,  and  diagrams. 
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II.  Theory 


Chapter  Overview 

The  purpose  of  this  chapter  is  to  review  the  theory  associated  with  this  thesis. 

The  chapter  is  divided  into  two  sections.  Section  one  covers  the  fundamental  physics  that 
make  the  Compton  Spectrometer  possible  and  some  of  the  fundamentals  of  radiation 
detection.  The  topics  include  photon  interactions  with  matter,  detector  response 
functions,  detector  resolution  and  efficiency,  and  attenuation.  Section  two  covers  the 
operation  of  the  Compton  Spectrometer  and  includes  system  configuration,  advantages 
and  disadvantages. 

Fundamentals  of  Nuclear  Physics  and  Radiation  Detection 

Photons  Interactions  with  Matter 

There  are  many  categories  of  photon  interactions  with  matter.  The  three  major 
interactions  are  photoelectric  absorption,  Compton  scattering  and  pair  production.  As 
shown  in  Figure  3,  which  effect  is  most  likely  to  occur  can  be  generalized  based  on  the 
energy  of  the  incoming  gamma  ray  and  the  atomic  number  (Z)  of  the  absorbing  material. 
The  line  separating  each  region  depicts  the  equal  probability  of  the  adjacent  effects 
occurring.  The  photoelectric  effect  (x)  is  dominant  at  low  gamma  ray  energies  and 
increases  with  higher  Z.  Compton  scattering  (a)  is  the  dominant  effect  for  low  Z 
absorbers  for  a  wide  range  on  energies.  Pair  production  (k)  requires  minimum  gamma 
ray  energy  of  1 .02  MeV  and  is  not  significant  for  energies  below  several  MeV.  This 
project  deals  with  lower-energy  photons,  therefore  pair  production  is  not  be  discussed 
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further.  Two  detector  materials  used  in  this  research,  cadmium  telluride  (CdTe)  and 
thallium  doped  sodium  iodide  (Nal(Tl)),  have  Zs  of  approximately  50  and  53 
respectively.  Figure  3  shows  that  for  these  materials  the  photoelectric  effect  and 
Compton  effect  are  comparably  efficient  for  gamma  rays  of  300  keV.  Over  the  majority 
of  the  energy  range  of  interest,  200-700  keV,  the  Compton  effect  is  dominant. 


Figure  3.  Relative  importance  of  the  three  major  types  of  gamma  ray  interactions 
based  on  photon  energy  and  atomic  number  (Knoll,  2000:52) 


Photoelectric  Effect 

Photoelectric  absorption  is  the  process  in  which  the  incoming  gamma  ray  is 
completely  absorbed  in  an  atom  resulting  in  an  electron  being  ejected.  The  photoelectron 
is  isotropic  and  has  the  same  energy  of  the  incident  gamma  ray,  less  the  binding  energy 
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of  the  ejected  electron.  If  there  is  sufficient  energy,  one  of  the  tightly  bound  K-shell 
electrons  can  be  liberated.  Characteristic  x-rays  may  then  be  emitted  as  outer  shell 
electrons  cascade  down  to  fill  lower  shell  vacancies.  No  simple  fonnula  for  the 
probability  of  photoelectric  absorption  exists  but  can  be  generalized  as  r  oc  Z4/ (h  v)3 , 
where  hv  is  the  gamma  ray  energy  (Turner,  1995:173). 


Compton  Scattering 

Compton  scatter  occurs  when  the  incident  gamma  ray  interacts  with  an  electron  in  the 
absorbing  material.  As  shown  in  Figure  4,  the  incident  gamma  ray  of  energy  h  v  is 
scattered  at  an  angle  0,  with  a  resulting  energy  of  h  v' ,  expressed  as: 

hv 

h  v  = - - - 

l  +  _^(l_cos*) 
m0c 

where  m0c2  is  the  rest  mass  energy  of  the  electron  (511  keV)  (Knoll,  2000:5 1). 


Figure  4.  Compton  scatter  diagram;  incident  photon  scattered  at  angle  0,  Compton 
electron  scattered  at  angle  <p  (Knoll,  2000:51) 
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The  energy  of  the  scattered  electron,  also  known  as  the  Compton  electron,  is  equal  to  the 
incident  gamma  ray  energy  minus  the  scattered  gamma  ray  energy  and  is  given  as 
(Turner,  1995:178): 


T 


=  hv- 


l-cos# 


m,c 
h  v 


+  1-COS# 


(2) 


The  Compton  electron  is  scattered  at  an  angle  (cp)  ranging  from  0-90  degrees  in  the 
opposite  direction  from  the  scattered  gamma  ray  in  order  to  conserve  momentum.  The 
relation  between  cp  and  9  is  given  as  (Turner,  1995:  178): 


cot 


0 


tan  cp  =  ■ 


1  + 


muc 

hv 


(3) 


The  Klein-Nishina  formula  calculates  the  angular  distribution  of  scattered  gamma  rays: 


da_  =  Zr 2 

dn  z,° 


l  +  o^l-cos#) 


V 


1  -  cos2  o) 
2 


1  + 


a2(l-cos61)2 
(l  +  cos  -COS61)) 


(4) 


where  a  =  hvjm^c  and  r0  is  the  classical  electron  radius  (Knoll,  2000:51).  Figure  5 

shows  the  Klein-Nishina  formula  plotted  for  varying  incident  gamma  ray  energies  from  1 
keV  to  5  MeV.  The  higher-energy  gamma  rays  clearly  tend  to  forward  scatter  while  very 
low-energy  gamma  rays  have  an  equal  probability  of  back  scatter. 
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Figure  5.  Polar  plot  of  Klein-Nishina  formula  representing  the  probability  of  the 
incident  gamma  ray,  of  initial  energy  shown,  scattering  at  an  angle  0 

Integrating  the  energy  dependent  Klein-Nishina  formula  over  all  solid  angles 
gives  the  total  Compton  cross  section.  The  fonnula  assumes  that  the  interaction  is  with 
free  electrons  and  does  not  consider  the  bound  state  of  the  electrons.  This  assumption 
leads  to  larger  Compton  cross  sections  for  lower-energy  gamma  rays.  For  CdTe,  the 
relative  error  of  the  calculated  value  to  experimental  data  range  from  1 1.5%  for  100  keV 
gamma  rays  to  1.1%  for  700  keV  gamma  rays  as  shown  in  The  integration  can  also  be 
conducted  over  a  portion  of  the  solid  angle  to  determine  the  probability  that  the  gamma 
ray  will  scatter  into  the  solid  angle  of  interest. 
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Table  1.  Calculations  are  shown  in  Appendix  C.  The  integration  can  also  be 
conducted  over  a  portion  of  the  solid  angle  to  detennine  the  probability  that  the  gamma 
ray  will  scatter  into  the  solid  angle  of  interest. 


Table  1.  Calculated  and  Experimental  Values  for  Compton  Cross-Section 
_ _ (Nowotny,  1998)  _ 


Incident  y 

Klein-Nishina 

XMuDat 

Relative  Error 

(keV) 

•  2\ 

(Centimeter  ) 

(%) 

100 

0.12360 

0.11090 

11.5 

200 

0.10200 

0.09727 

4.9 

300 

0.08868 

0.08628 

2.8 

400 

0.07946 

0.07782 

2.1 

500 

0.07256 

0.07150 

1.5 

600 

0.06712 

0.06624 

1.3 

700 

0.06268 

0.06200 

1.1 

Attenuation  and  Probability  of  Interaction 

The  attenuation  of  gamma  rays  can  be  expressed  as  a  function  of  the  linear 
attenuation  coefficient  ( ju  )  and  the  thickness  of  absorber  material  (t).  The  sum 
probabilities  of  the  individual  interaction  processes  gives  jU  =  t  +  <j  +  ic  .  Given  an  initial 
number  of  gamma  rays  emitted  ( I0 )  and  a  transmitted  number  of  gamma  rays  (I)  the 
attenuation  is  expressed  as  (Tsoulfanidis,  1995:  159) 

I=I0e-r  (5) 

The  probability  of  an  interaction  occurring  between  the  depths  x1  and  x2  as  shown  in 
Figure  6  can  then  be  expressed  as 
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Figure  6.  Probability  of  gamma  ray  interaction  between  depths  xl  and  x2  expressed 
as  the  difference  of  attenuation  at  xl  and  x2 

Detector  Response  to  Gamma  Radiation 

The  response  of  a  detector  can  be  described  by  the  probability  that  electrons  and 
photons,  “produced  in  interactions  of  the  original  gamma  rays,”  escape  the  detector 
volume  before  losing  all  of  their  energy  (Knoll,  2000:  3 12).  The  electrons  and  photons 
are  created  from  photoelectric  and  Compton  interactions  and  include  the  photoelectrons, 
Compton  electrons,  Compton  scattered  gamma  rays  and  bremsstrahlung  radiation.  The 
detector  response  and  can  be  generalized  in  three  size  categories  of  small,  intennediate 
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and  very  large  (Knoll,  2000:3 12-317).  The  size  refers  to  the  detector’s  dimensions 
compared  to  the  mean  free  path  of  subsequent  electrons  and  photons.  In  a  very  large 
detector  the  mean  free  path  is  small  compared  to  the  detector  dimensions  and  all  energy 
is  converted  to  information  carriers.  All  interactions  therefore  correspond  to  counts 
under  the  “full-energy  peak.”  In  small  detectors  the  dimensions  are  small  compared  to 
the  mean  free  path  of  the  electrons  and  photons  and  energy  is  lost  when  they  escape. 
Therefore  fewer  information  carriers  are  created  and  counts  appear  below  the  full-energy 
peak.  In  an  intermediate  size  detector  the  probability  of  subsequent  photons  interacting 
before  escape  increases  resulting  in  more  information  carriers  created.  This  increases  the 
number  of  counts  closer  to  the  full-energy  peak. 

The  Amptek  CdTe  detector  used  has  dimensions  of  1.0  x  5.0  x  5.0  mm3  and 
generally  displays  small  detector  characteristics.  In  this  research  the  gamma  ray 
interaction  of  interest  is  a  Compton  event  resulting  in  a  scattered  gamma  ray  of 
approximately  500  keV  escaping  the  detector  and  a  Compton  electron  of  approximately 
100  keV  being  collected  in  the  detector.  In  CdTe  the  mean  free  path  of  a  500  keV 
gamma  ray  is  1.8  centimeters  and  the  range  of  a  100  keV  electron  is  0.03  mm  (Williams, 
2004:  30).  The  relative  size  of  the  mean  free  paths  to  the  detector  dimensions  shows  that 
the  scattered  gamma  ray  has  a  high  probability  of  escaping  the  detector  volume  while  the 
Compton  electron  does  not.  Materials  research  continues  to  improve  CdTe  growth 
techniques  but  the  technology  for  making  larger  CdTe  detectors  with  good  energy 
resolution  does  not  yet  exist.  The  current  size  limitation  is  one  of  the  reasons  this 
research  is  attractive  and  possible. 
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As  shown  in  Figure  7,  the  incident  gamma  ray,  y  (1),  Compton  scatters  in  the 
detector  with  the  scattered  gamma  escaping  while  the  Compton  electron  is  completely 
absorbed.  Based  on  the  Klein-Nishina  formula  and  the  Compton  relation  shown  in 
Equations  2  and  3,  a  distribution  of  Compton  scattered  electron  energies,  called  the 
“Compton  continuum,”  is  characteristic  of  small  detectors.  The  maximum  energy 
deposited  by  a  single  Compton  scatter  corresponds  to  the  gamma  ray  backscattering  with 
0=7i  and  is  given  as  (Knoll,  2000:3 1 1): 


E,  = 


hv 


1  +  2A  W  m0c 


(7) 


Figure  7.  Various  gamma  ray  interactions  affecting  detector  response 

This  is  known  as  the  “Compton  edge”  and  is  shown  in  Figure  8.  In  intermediate  size 
detectors,  it  is  possible  for  multiple  Compton  scatters  to  occur,  as  illustrated  with  y  (2),  in 
Figure  7.  This  response  leads  to  counts  between  the  Compton  edge  and  the  full-energy 
peak. 
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Figure  8.  Calculated  Compton  continuum  and  edge  for  a  662  keV  photon 
in  CdTe  (full-energy  peak  not  shown) 

Although  the  interaction  of  secondary  gamma  radiation  as  described  above  is 
unlikely  in  a  small  detector  it  can  be  more  prevalent  based  on  the  source-detector 
geometry.  Figure  9  shows  two  source-detector  geometries.  The  dashed  arrows  represent 
a  typically  forward  scattered  gamma  ray  emitted  from  a  Compton  scatter  event  at  a 
specific  angle  and  the  mean  free  path  of  the  scattered  gamma  ray.  Because  of  the  larger 
angle  created  by  the  position  of  Source  A,  scattered  gamma  ray  y’  ia  must  travel  farther 
through  the  detector  material  than  y’2A,  y’m,  or  y’2B  to  escape  the  detector  volume. 
Therefore  the  Source  A-detector  geometry  has  a  higher  probability  of  producing  multiple 
Compton  scatter  events  than  the  Source  B-detector  geometry. 
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Detector 


Figure  9.  Example  of  source-detector  geometry  effects  on  detector  response 
Detector  Response  Complications 

The  expected  detector  response  can  also  be  distorted  because  of  secondary 
electrons  or  bremsstrahlung  radiation  escaping  the  detector  volume  before  being  absorbed 
as  shown  by  y  (3)  and  y  (4),  in  Figure  7.  In  small  size  detectors  the  possibility  of  the 
secondary  electrons  leaking  from  the  detector  surface  increases.  Electron  loss  also 
increases  for  higher-energy  incident  radiation  because  the  secondary  electrons  generally 
have  higher  energy  and  therefore  longer  range  in  the  detector.  This  leakage  results  in  the 
Compton  spectrum  to  be  altered  favoring  lower  energies.  Bremsstrahlung  radiation  is 
electromagnetic  radiation  created  as  the  secondary  electrons  lose  energy  in  the  detector. 
Turner  expresses  the  ratio  of  bremsstrahlung  (radiative)  to  collisional  stopping  powers  for 
an  electron  of  energy  E  (MeV),  in  a  material  of  atomic  number  Z  as 
Radiative/ Collisional  =  ZE/ 800  (Turner,  1995:  143).  In  CdTe  the  stopping  power  ratio 
for  662  and  100  keV  electrons  is  .04  and  .006  respectively.  It  is  clear  that  bremsstrahlung 
is  not  a  significant  factor  in  the  response  of  the  detector. 
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Detector  Resolution  and  Efficiency 


Detector  resolution  describes  how  well  a  detector  can  distinguish  an  energy  peak 
of  a  given  incident  radiation  and  is  quantified  using  the  full-width  at  half-maximum 
(FWHM)  divided  by  the  peak  centroid  energy  H0.  Similarly  the  resolution  can  also  be 
described  using  the  full-width  at  tenth-maximum  (FWTM)  divided  by  the  peak  centroid 
energy  Ho.  These  standards  are  based  on  a  Gaussian  shaped  peak.  Figure  10  shows  an 
example  spectrum  for  88  keV  photons  and  the  values  used  to  calculate  resolution. 


Figure  10.  Example  energy  spectrum  for  full- width  at  half-maximum 
and  full-width  at  tenth-maximum  resolution  calculations 
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The  calculated  values  from  Figure  10  for  Rfwhm  and  Rfwtm  are  1.3%  and  5.0%, 
respectively.  The  large  difference  in  FWHM  and  FWTM  of  this  example  is  characteristic 
for  some  semiconductor  detectors.  The  asymmetric  peak  is  caused  by  incomplete 
information  carrier  collection  or  secondary  electron  and  bremsstrahlung  escape  from  the 
detector  (Knoll,  2000:  434-435).  Consequently  careful  consideration  of  the  detector 
system  and  application  must  be  taken  when  choosing  a  resolution  calculation  method. 

Knoll  describes  three  factors  that  contribute  to  the  FWHM  (Wt)  in  semiconductor 
detectors;  charge  carrier  statistics  (WD),  charge  carrier  collection  (Wx)  and  electronic 
noise  (WE)  (Knoll,  2000:  417-419).  The  sum  of  the  squared  factors  gives  the  squared 
total  FWHM  as  follows 


W?  =W£  +  Wt  +  Wz. 


(9) 


The  first  factor  is  due  to  the  statistical  fluctuation  in  charge  carriers  created  and  is 
calculated  as 


Wo  ={2.35)2FsE,  (TO) 

where  F  is  the  Fano  factor,  s  is  the  energy  needed  to  create  an  electron-hole  pair,  and 
E  is  incident  gamma  ray  energy.  The  Fano  factor  is  used  to  correct  for  observed 
deviations  from  Poisson  statistics  and  ranges  from  unity  for  scintillation  detectors  to  less 
than  0. 1  for  small  germanium  detectors  (Knoll,  2000:  1 15-116).  The  Fano  factor  for 
CdTe  is  0.15  and  s  is  4.5  eV  (Takahashi,  2000:  2). 
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The  second  contribution,  Wx,  is  due  to  incomplete  charge  carrier  collection  which 
occurs  when  charge  carriers  are  trapped  or  recombine  before  being  collected.  The  effects 
of  incomplete  charge  collection  increase  with  larger  detector  volumes  and  lower  electric 
fields.  The  contribution  can  be  found  experimentally  by  conducting  FWHM 
measurements  while  varying  the  electric  field.  The  assumption  is  that  the  effects  of 
incomplete  charge  collection  are  insignificant  when  an  infinitely  large  electric  field  is 
applied  (Knoll,  2000:  417). 

The  last  factor,  We,  is  due  to  electronic  components  following  the  detector.  It  is 
measured  directly  by  finding  the  FWHM  of  the  spectrum  produced  by  a  precision  pulse 
generator  connected  to  the  preamplifier.  A  parallel  test  point  is  often  provided  for  this 
measurement  (Knoll,  2000:  418). 

The  ability  to  distinguish  two  closely  spaced  peaks  is  essential  to  determining  the 
isotopic  ratio  of  a  plutonium  sample  and  is  quantified  using  the  figure  of  merit,  M,  shown 
in  Figure  1 1 .  The  figure  of  merit  is  given  as 


M  = 


X 

W„  +  W,  ’ 

a  b 


(ID 


where  Wa  and  Wb  are  the  FWHM  of  each  peak  and  X  is  the  difference  between  the 

peaks’  centroids.  Generally,  a  minimum  value  of  one  for  M  is  needed  to  distinguish  two 
peaks,  and  larger  values  are  preferred. 
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Figure  11.  Figure  of  Merit  describes  the  ability  of  a  detector  to  distinguish  two 

energy  peaks  (Knoll,  2000:  680) 

Detector  efficiency  describes  the  ratio  of  pulses  created  by  the  detector  to  a 
number  of  radiation  quanta  emitted  from  a  source  and  is  critical  in  low  count  rate 
applications  such  as  treaty  monitoring.  Absolute  efficiency  considers  the  total  number  of 
quanta  emitted  from  a  source  while  intrinsic  efficiency  considers  only  the  radiation 
quanta  incident  on  the  detector  volume.  For  an  isotropic  source  they  are  related  as 

=  £'a/„  (47r/0) ,  where  Q  is  the  solid  angle  from  the  source  subtended  by  the  detector. 

The  solid  angle  depends  on  the  source-detector  geometry.  For  parallel  planar  source  and 
detector  geometry,  see  Figure  12,  the  solid  angle  is  calculated  as 

\\(dAjAw2)fiAd{n-rlr) 


Q  = 


A, 


(12) 


where  As  and  Ad  are  the  source  and  detector  areas  respectively,  dAs  and  dAd  are  the 
differential  source  and  detector  areas  respectively,  r  is  the  vector  from  dAd  to  dAs  with 
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magnitude  r ,  and  n  is  a  unit  vector  normal  to  detector  surface  (Tsoulfanidis,  1995:  268- 
269). 


Figure  12.  Parallel  planar  detector-source  geometry  for  calculation  of  solid  angle 
subtended  by  detector  from  source  (Tsoulfanidis, 1995:  269) 

Tsoulfanidis  presents  a  solid  angle  approximation  for  disk  source  and  rectangular 

detector  geometries: 

Q  =  C°^°2  +o22)+^-(5^4  +3co*  +  3a>2')+‘^Y2if°\  +<^22) 

-—i//6  +—co2a>l  y/4{<x>l  +a>2  )-—i//2(9ci>\  +9 cot  +I0a>2a>2)  (13) 

64  12  16  v  7  32  v  ’ 

- (0,  ( 0  [(')  +  (0  ) - 1 (')  +  ) I 

16  i  -v  1  2)  i6 v  i  2 

where  a  and  b  are  the  rectangular  detector  dimensions,  d  is  the  source  distance  above 
the  detector,  Rs  is  the  source  radius,  o\  =  a/ d ,  co2  =  bl d  ,  and  i//  =  Rs  / d  .  This 

approximation  is  valid  for  geometries  where  cox ,  <x>2 ,  and  i//  are  less  than  one 
(Tsoulfanidis,  1995:276-275). 
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Gotoh  and  Yagi  used  inverse  trigonometric  functions  to  determine  the  solid  angle 


subtended  by  a  rectangle  at  an  arbitrary  point  [xp ,  y  ,  z  ) : 
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where  the  rectangular  detector  is  enclosed  by  the  four  lines  x  =  x1 ,  x  =  x2,  y  =  y, ,  and 
y  =  y2  (Gotoh,  1971:  485). 


Signal-to-Noise  Ratio  and  Minimum  Detectable  Activity 

“Noise  is  any  undesired  fluctuation  that  appears  superimposed  on  a  signal  source” 
(Knoll,  2000:  629).  Noise  affects  the  resolution  of  all  detector  systems  but  becomes 
critical  when  the  noise  variance  is  of  the  same  or  larger  order  of  magnitude  of  the  source 
variance.  Generally,  noise  affects  scintillation  detectors  less  than  semiconductor 
detectors.  Where,  in  the  signal  chain  the  noise  occurs,  significantly  affects  the  final 
magnitude  of  the  noise.  Noise  at  the  beginning  is  amplified  the  same  as  the  true  signal 
therefore  electronic  noise  often  focuses  on  the  preamplifier  and  its  input  stage.  The 
signal  to  noise  ratio  can  be  quantified  as  S  /  N  =  Ns  /NB  ,  where  Ns  =  NT  -  NB  is  the 
counts  attributed  to  the  source,  NT  is  the  total  counts  recorded,  and  N B  is  the  counts 
attributed  to  background  noise. 

The  minimum  detectable  amount  (MDA)  represents  the  minimum  measurable 
activity  based  on  the  detection  system  to  provide  a  desired  level  of  certainty  that 
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radioactive  materials  are  truly  present.  Equation  15  is  often  called  the  “Curie  Equation” 
and  gives  the  minimum  number  of  counts  N  D  to  ensure  less  than  5%  false-negative 

counts  and  5%  false-positive  given  ctNb  =  ^NB  (Knoll,  2000:  94-96). 


Nd  =4.653(7^  +2.706  (15) 

A  minimum  detectable  activity  ( a  )  can  then  be  calculated  using  the  radiation  yield  per 
disintegration  ( /  ),  the  absolute  detection  efficiency  ( e  )  and  the  counting  time  ( T )  and 
is  given  as  (Knoll,  2000:  96) 


a  = 


No_ 

M 


(16) 


The  Compton  Spectrometer 

The  Compton  spectrometer  consists  of  two  detectors  operating  in  a  coincidence 
mode  as  shown  in  Figure  13.  The  incident  gamma  rays  enter  the  first  detector  where  they 
interact  according  to  the  linear  attenuation  coefficient  as  discussed  above.  Compton 
scatter  events  produce  Compton  electrons  that  are  collected  in  the  detector,  referred  as  the 
electron  energy  analyzer  (EEA).  A  fraction  of  the  Compton  scattered  gamma  rays  will 
reach  the  second  detector  referred  to  as  the  photon  energy  analyzer  (PEA).  Because  of 
the  short  distances  traveled,  the  electron  and  photon  pulses  created  are  essentially 
coincident  events.  Associated  electronics  detennine  which  events  are  true  coincidence 
and  the  energy  spectra  of  both  detectors  are  collected.  The  EEA  spectrum  will  consist  of 
Compton  scattered  electron  energies  given  by  Equation  2  with  9  ranging  from  9mm  to 
#max  •  Similarly,  the  PEA  spectrum  consists  of  the  scattered  gamma  ray  energies  given  by 
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Equation  1  with  0  ranging  from  0mm  to  6mm .  Some  variation  in  6min  and  #max  will  occur 
due  to  the  point  of  interaction  in  the  EEA. 


Figure  13.  The  Compton  spectrometer  consisting  of  the  electron  energy  analyzer 
(EEA)  and  the  photon  energy  analyzer  (PEA) 

The  spectra  collected  in  both  the  EEA  and  PEA  are  affected  by  the  characteristic 
efficiency  and  resolution  typical  of  those  detector  types.  The  coincidence  requirement  of 
the  system  significantly  reduces  background  counts  and  therefore  should  increase  the  S/N 
ratio.  Based  on  the  expected  incident  gamma  ray  energies,  the  EEA  is  selected  to  provide 
the  best  resolution  while  maximizing  the  number  of  Compton  scatter  events.  As 
previously  discussed,  the  Compton  cross  section  dominates  for  a  wide  range  of  energies 
especially  in  low  Z  materials.  The  PEA  is  selected  for  high  efficiency,  while  resolution  is 
less  of  a  concern. 
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There  are  several  advantages  and  disadvantages  associated  with  the  Compton 
spectrometer.  The  advantages  include  an  improved  S/N  ratio  from  reducing  background 
counts  and  better  resolution  from  collecting  the  forward  scattered  and  lower-energy 
Compton  electron.  As  previously  mentioned,  the  background  radiation  is  expected  to  be 
significantly  reduced  because  of  the  coincidence  requirements.  Similarly,  the  spectral 
artifacts  such  as  backscatter  peaks  and  x-ray  escape  peaks  are  also  reduced  because  of  the 
coincidence  requirements.  The  improved  S/N  ratio  is  expected  to  significantly  reduce  the 
noise  contribution  to  the  FWHM.  The  collection  of  the  Compton  electron  versus  the 
photoelectron  also  reduces  the  energy-dependent  contributions  to  the  FWHM. 

The  Wx  and  Wd  contributions  to  the  FWHM  increase  with  the  amount  of  energy 
deposited  in  the  detector,  therefore  the  lower-energy  Compton  electrons  are  collected  at  a 
higher  resolution.  The  Compton  electrons  generated  from  incident  gamma  rays  are 
always  lower  in  energy  than  the  photoelectrons  generated  from  the  same  gamma  rays. 

For  small  gamma  ray  scatter  angles  the  Compton  electron  energy  can  reach  an  order  of 
magnitude  less  than  the  photoelectrons.  The  higher  resolution  Compton  electron 
spectrum  collected  in  the  EEA  can  then  be  used  to  produce  an  incident  gamma  ray 
spectrum  based  on  the  non-linear  Compton  relation  as  demonstrated  in  Figure  14.  In  the 
example  spectra,  the  100-keV  Compton  electron  peak  has  a  resolution  of  0.96%,  which  is 
reasonably  attainable  for  this  energy.  The  calculated  662-keV  incident  gamma  ray  peak 
has  a  resolution  of  0.52%.  This  resolution  is  difficult  to  achieve  without  cryogenics. 
Calculations  and  further  analysis  of  the  Compton  relations  are  found  in  Appendix  D. 

Also  an  advantage  of  collecting  the  Compton  electrons  is  their  predominance  of 
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scattering  in  the  forward  direction  from  0-90  degrees  (Turner,  1995:  178)  resulting  in  less 
probability  of  the  electrons  escaping  small  detector  volumes  before  complete  energy 
collection. 
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Figure  14.  Demonstration  of  the  partial-energy  Compton  electron  spectrum 
collected  in  the  EEA  converted  to  the  calculated  full-energy  gamma  ray  spectrum 
based  on  Compton  relation;  Improved  resolution  from  0.96%  to  0.52% 

The  main  disadvantage  of  the  system  is  efficiency.  The  total  Compton  cross 
section  in  CdTe  is  predominant  over  the  300  to  700  keV  gamma  ray  energy  range.  In  the 
Compton  spectrometer  the  solid  angle  is  subtended  by  the  PEA  is  only  a  small  portion  of 
the  47i  steradians  used  to  detennine  the  total  Compton  cross  section.  Therefore,  only  a 
small  portion  of  the  total  Compton  cross  section  is  utilized.  For  example,  the  total 
Compton  cross  section  for  a  662-keV  photon  in  CdTe  is  6.4x10'  cm7gram  while  the 
Compton  cross  section  utilized  in  the  Compton  spectrometer  is  1 .9x10  4  cm2/gram.  This 
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is  nearly  two  orders  of  magnitude  less  than  the  photoelectric  cross  section  of  1.2x1  O'2 
cm  /gram  for  the  same  photon  energy.  The  efficiency  problem  is  further  compounded  by 
strict  source  collimation  requirements.  To  maintain  a  reasonable  resolution  in  the 
Compton  electron  spectrum,  the  source  must  be  collimated  to  the  EEA  to  minimize  the 
scatter  angle  deviations.  Collimation  and  alignment  considerations  of  the  system  play  a 
significant  role  in  the  effectiveness  of  the  system. 
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III.  Equipment 


Chapter  Overview 

The  purpose  of  this  chapter  is  to  discribe  the  equipment  used  in  this  thesis.  The 
chapter  is  divided  into  five  sections  covering  the  complete  system  configuration, 
detectors,  electronics,  collimators  and  radioactive  sources.  First,  an  overview  of  the 
system  configuration  is  discussed.  Then  the  Amptek  XR-IOOT-Cadmium  Telluride 
(CdTe)  and  Bicron  Model  3M3/3  thallium-doped  sodium  iodide  (Nal(Tl))  detectors  and 
their  associated  pulse-processing  electronics  are  reviewed.  Next,  the  purpose  and  general 
characteristics  of  the  nuclear  instrument  module  (NIM)  pulse-processing  electronics  used 
for  coincidence  determination  and  spectrum  collection  are  covered,  followed  by  an 
examination  of  the  materials  and  dimensions  of  the  linear  and  conical  collimators. 

Finally,  the  radioactive  sources  used  are  outlined. 

Compton  Spectrometer  System 

The  Compton  spectrometer  system  used  in  this  thesis  is  based  on  the  design  used 
by  Capt.  Williams  shown  in  Figure  15  (Williams,  2004:  44-46).  The  gamma  rays 
originate  at  position  1  from  a  typically  isotropic  sample.  As  they  move  through  position 
2,  the  gamma  rays  are  collimated  to  the  EEA  (CdTe)  at  position  3  and  more  importantly, 
they  are  aligned  normal  and  centered  on  the  conical  collimator  at  position  4.  The  gamma 
rays  interact  in  the  EEA  and  some  scatter  out  through  the  conical  collimator  to  the  PEA 
(Nal(Tl))  at  position  5.  When  coincidence  pulses  occur  in  the  pulse-processing 
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electronics  the  energy  data  of  the  EEA  and  PEA  are  sent  to  the  spectrum-collection 
electronics. 


1  ■ 


Figure  15.  Compton  spectrometer  configuration;  1.  Isotropic  sample, 

2.  Initial  linear  collimator,  3.  EEA  (CdTe  detector),  4.  Conical  collimator, 

5.  PEA  (Nal(Tl)  detector) 

Detectors 

Electron  Energy  Analyzer 

The  EEA  consists  of  the  Amptek  Model  XR-lOOT-CdTe  detector  and  the  Model 
PX2T  power  supply  and  amplifier  which  are  pictured  and  tabulated  in  Appendix  E.  The 
5x5x1  mm  CdTe  crystal  is  mounted  on  a  thennoelectric  cooler  as  shown  in  Figure  16. 
The  cooler  maintains  the  detector  crystal,  the  field-effect  transistor  (FET),  and  feedback 
components  at  approximately  -30  degrees  Celsius,  which  can  be  monitored  through 
connections  at  the  PX2T  (Amptek,  2002:  20). 
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Figure  16.  CdTe  hermetic  detector  housing  (Amptek,  2002:  20) 

The  detector  housing  is  designed  for  gamma  rays  to  enter  through  a  10  mil 
beryllium  window  incident  on  the  5x5  mm  CdTe  crystal  face.  In  the  Compton 
spectrometer  configuration,  the  gamma  rays  enter  through  the  10-mil  thick,  nickel  wall  of 
the  detector  housing  to  strike  the  1x5  mm  edge  of  the  CdTe  crystal.  The  difference  in 
gamma  ray  attenuations  for  nickel  and  beryllium  has  minimal  effect  on  the  system  as 
shown  in  Appendix  F.  The  output  of  the  FET  is  sent  to  the  Amptek  A250  charge 
sensitive  preamplifier.  The  negative  preamplifier  signal  is  sent  to  the  PX2T.  The  PX2T 
provides  the  +400  volt  detector  bias,  the  ±8  volt  preamplifier  power,  the  +8  volt 
temperature  monitor  power  and  the  0-3  volt  cooler  power.  The  PX2T  also  includes  a 
triangular  pulse-shape  amplifier,  output  shown  in  Figure  17,  and  a  rise-time  discriminator 
(RTD).  The  RTD  reduces  the  characteristic  “hole  tailing”  of  the  energy  spectrum  created 
by  inefficient  hole  collection  which  degrades  energy  resolution  as  shown  in  Figure  18. 
Essentially  the  RTD  discards  pulses  that  originate  from  interactions  farther  into  the 
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detector  from  where  the  holes  must  travel  a  large  portion  of  the  detector  width  before 
being  collected. 
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Figure  17.  PX2T  Shaping  Amplifier  output  with  RTD  inactive  (Amptek,  2002:  21) 

This  longer  travel  time  increases  the  probability  of  recombination  resulting  in  lower- 
energy  pulses  recorded.  The  improvement  in  resolution  with  the  RTD  active  is  at  the  cost 
of  efficiency  since  the  active  volume  of  the  detector  is  decreased.  The  RTD  is  activated 
at  the  PX2T,  and  the  sensitivity  can  be  adjusted  internally  at  R22.  To  evaluate  the 
decrease  in  efficiency  the  PX2T  provides  a  rear  panel  input  count  rate  (ICR).  A  lower 
level  discriminator  is  set  just  above  the  noise  level  creating  a  short  output  pulse  for  all 
events  regardless  of  energy-collection  time.  This  count  rate  can  then  be  compared  with 
the  count  rate  from  the  RTD  to  determine  a  relative  efficiency.  In  the  characterization  of 
the  detector  a  Canberra  Model  3 125  high  voltage  power  supply  (HVPS)  is  used  to 
provide  the  variable  bias  voltage  for  the  CdTe  detector. 
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Figure  18.  Energy  spectrum  with  and  without  RTD  active  demonstrating  the 
reduction  in  hole  tailing  (Amptek,  2002:  21) 

Photon  Energy  Analyzer 

The  PEA  consists  of  the  Bicron  Model  3M3/3  thallium-doped  sodium  iodide 
(Nal(Tl))  scintillation  detector  and  pulse-processing  electronics  pictured  and  specified  in 
Appendix  G.  The  Nal(Tl)  crystal  has  a  3  inch  diameter  and  3  inch  depth  resulting  in  a 
2 1 .2  inch3  volume.  “Nal(Tl)  produces  the  highest  signal  in  a  photomultiplier  tube  (PMT) 
per  amount  of  radiation  absorbed  in  the  crystal  of  all  presently  known  scintillators” 
(Bicron,  2002:  1).  The  Bicron  PA- 14  preamplifier  connects  directly  to  the  14-pin  output 
of  the  PMT  and  provides  gain  balance  (G)  and  focus  (F)  potentiometers  for  fine 
adjustments  in  pulse  height  and  resolution.  The  detector  is  operated  at  a  bias  voltage  of 
approximately  +900  volts,  and  is  provided  by  a  computer-controlled  Canberra  9645 
FIVPS.  The  -24-volt  power  required  for  the  PA- 14  is  a  standard  output  on  ORTEC 
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amplifiers.  Both  the  ORTEC  Model  672  Spectroscopy  Amplifier  and  ORTEC  579  Fast 
Filter  Amplifier  were  used  to  process  the  Nal(Tl)  energy  pulses.  The  672  provided  more 
shaping  options  and  better  energy  resolution  with  less  stable  timing  characteristics 
(ORTEC,  2004:  3.46-3.49).  On  the  other  hand,  the  579  provides  a  fast  rise  time  and 
better  timing  characteristics  (ORTEC, 2004:  3.36-3.39)  while  sacrificing  energy 
resolution.  Bicron  quoted  specification  for  the  Nal(Tl)  detectors  is  7.0  ±0.2%  resolution, 
resulting  in  a  FWHM  of  less  than  50  keV  for  662  keV  gamma  rays. 

Coincidence  Electronics 

Two  methods  for  detennining  the  coincidence  of  the  EEA  and  PEA  pulses  were 
used;  the  pulse  start-stop  technique  using  a  time-to-amplitude  converter  (TAC)  and  the 
pulse  overlap  technique  using  a  slow  coincidence  unit.  The  equipment  used  for  both 
methods  is  specified  in  Appendix  H. 

In  both  methods  the  first  step  is  to  discriminate  the  timing  of  a  pulse’s  arrival. 

This  process  is  complicated  by  amplitude  and  rise-time  walk  and  jitter  as  shown  in  Figure 
19.  A  timing  pulse  is  generated  when  a  trigger  level  is  crossed,  but  the  timing  of  two  true 
coincident  events  can  be  different  based  on  the  pulse  shapes.  Most  ORTEC  electronics 
use  a  trailing-edge  constant-fraction  method  to  minimize  the  effects  of  amplitude  walk 
but  are  still  susceptible  to  variations  in  pulse  shapes. 
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Figure  19.  Amplitude  walk,  rise-time  walk  and  jitter  effects  on  pulse  discrimination 

(Knoll,  2002:  659-661). 

The  ORTEC  Model  551  Timing  Single-Channel  Analyzer  (T-SCA)  and  Model 
552  Pulse-Shape  Analyzer  (PSA)/T-SCA  were  used  to  derive  the  timing  signals.  Both 
modules  use  ORTEC ’s  “patented  trailing-edge  constant- fraction  (CF)  timing  technique” 
(ORTEC,  2004:  4.9)  with  optional  rear  panel  connectors  providing  standard  leading  edge 
timing.  The  main  difference  between  the  modules  is  the  55 1  ’s  CF  is  set  at  50%  while  the 
552  provides  an  A  and  B  channel  output  both  with  separate  CF  settings.  This  allows 
information  about  the  pulse-shape  to  be  obtained.  Each  module  can  provide  positive  or 
negative  NIM-standard  output  pulses  (ORTEC,  2004:  4.7-4. 1 1). 
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While  using  the  TAC  method,  the  fast-negative  timing  pulses  from  the  PSA/T- 
SCAs  are  sent  directly  to  the  ORTEC  Model  566  TAC.  The  TAC  provides  a  square 
output  pulse  with  amplitude  proportional  to  the  time  between  the  start  and  stop  pulses. 

To  limit  dead-time  associated  with  receiving  a  second  start  pulse  before  a  stop  pulse,  the 
higher-count-rate  output  of  the  Nal(Tl)  is  connected  to  the  stop  gate.  The  output 
spectrum  from  the  TAC  should  show  a  “prompt  coincidence  peak”  associated  with  true 
coincidence  on  top  of  a  chance  coincidence  continuum  (Knoll,  2002:  666-667).  A  SCA 
is  then  used  to  select  the  prompt  coincidence  peak. 

Operating  in  the  pulse  overlap  method,  square  logic  pulses  are  required  at  the 
input  of  the  ORTEC  Model  409  Linear  Gate  and  Slow  Coincidence  (LGSC)  unit.  The 
narrow  T-SCA  outputs  are  sent  to  ORTEC  Model  416A  Gate  and  Delay  Generators 
(GDG)  to  adjust  the  delay,  width,  and  amplitude  of  the  gating  pulses  (ORTEC,  2004: 

1 1.8).  The  LGSC  uses  the  simple  AND  gate  producing  a  logic  gate  pulse  anytime  the 
coincidence  inputs  overlap.  The  coincidence  alignment  and  resolving  times  are 
determined  by  the  delay  and  pulse  width,  respectively,  from  the  adjustable  GDGs. 

The  logic  gate  produced  from  the  LGSC  or  the  TAC-SCA  combination  is  then 
sent  to  the  Canberra  Model  9633  Analog-to-Digital  Converters  (ADC)  parallel  to  the 
energy  pulses  of  the  EEA  and  PEA.  The  ADC  processes  only  those  energy  pulses  with 
coincident  gate  pulses,  assigning  each  pulse  to  an  energy  bin  as  established  with  the 
Canberra  Genie  2000  software.  The  data  is  then  transferred  through  the  Canberra  Model 
5556A  Acquisition  Interface  Module  (AIM)  to  a  PC  running  the  Genie  2000  software.  A 
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single  AIM  is  capable  of  connecting  multiple  ADCs  to  a  single  PC  and  the  system  was 
established  to  collect  both  the  EEA  and  PEA  energy  spectra  simultaneously. 

Collimators 

The  Compton  spectrometer  requires  an  initial  linear  collimator  from  the  source  to 
the  EEA  and  a  conical  collimator  from  the  EEA  to  the  PEA  as  shown  in  Figure  15.  Two 
different  initial  collimators,  shown  in  Appendix  I,  were  used;  one  with  a  circular  cross 
section  of  radius  1 .5  mm  and  the  other  with  a  hexagonal  cross  section  with  internal  radius 
of  1.5  mm.  The  circular  collimator  is  made  of  lead  and  was  drilled  by  the  AFIT  model 
shop.  A  source  holder  was  centered  on  the  hole  and  attached  to  facilitate  consistent 
system  setup.  The  hexagonal  collimator  is  the  original  collimator  used  by  Capt.  Williams 
(Williams,  2004:  45).  The  conical  collimator  was  constructed  of  AIM70,  a  lead-bismuth 
alloy  with  a  melting  point  of  70°  Celsius,  by  Capt  Williams  (Williams,  2004:  122-123) 
and  its  dimensions  and  diagrams  are  located  in  Appendix  J.  The  angle  of  the  cone  is  0.53 
radians,  which  corresponds  to  a  100-keV  Compton-scattered  electron  from  an  incident 
662  keV  gamma  ray.  The  inside  cone  is  parallel  to  the  outside  cone  and  by  adjusting  its 
height  relative  to  the  outside  cone  changes  the  collimator  gap  size.  Spacers  are  used  to 
select  different  gap  sizes. 

Radioactive  Sources 

Many  radioactive  sources  were  used  throughout  this  thesis  and  a  tabulated  list  is 

57  137  152 

located  in  Appendix  K.  "  Co,  Cs,  Eu  and  multinuclide  sources  were  all  used  for 

22 

energy  calibration  and  detector  characterization.  The  “  Na  5 1  IkeV  annihilation  gamma 
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rays  were  used  for  aligning  the  EEA  and  PEA  pulses  for  coincidence.  137Cs  was  the 
primary  isotope  used  to  evaluate  the  Compton  spectrometer.  All  three  available  137Cs 


sources  were  combined  to  provide  an  activity  of  approximately  25  pCi.  The  adjusted 
activity  of  the  combined  source  ignored  self-absorption.  Adjusted  activity  calculations 
are  demonstrated  in  Appendix  L. 
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IV.  Methodology 


Chapter  Overview 

This  chapter  outlines  five  key  actions  taken  during  this  thesis  research  and  their 
purpose:  1.  test  and  characterize  the  CdTe  and  Nal(Tl)  detectors  individually,  2. 
construct  system  housing  to  facilitate  a  reproducible  geometry  with  improved  component 
alignment,  3.  establish  EEA  and  PEA  timing  pulse  alignment  and  assess  both  coincidence 
methods  described  above,  4.  collect  Compton  spectrometer  data  for  system  analysis  and 
5 .  develop  and  assess  a  computer  simulated  Compton  spectrometer  to  further  understand 
the  geometric  relations  and  interaction  probabilities  of  the  system. 

Detector  Testing  and  Characterization 

To  best  evaluate  the  Compton  spectrometer’s  performance,  the  individual 
detectors’  resolution  and  efficiency  were  detennined  first.  First,  the  CdTe  resolution  and 
efficiency  were  determined  using  the  geometry  shown  in  Figure  20.  The  system  was 
evaluated  with  the  rise  time  discriminator  (RTD)  on  and  off,  with  l37Cs  and  57Co  in 
positions  one  and  two.  Sample  57Co  pulse-height  spectra  for  each  geometry  are  shown  in 
Figure  21  and  Figure  22. 
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Position  2 


Figure  20.  CdTe  resolution  and  efficiency  geometry 


Figure  21.  Sample  57Co  spectrum  using  CdTe  detector 
(RTD  on,  Position  1) 


Figure  22.  Sample  57Co  spectrum  using  CdTe  detector 
(RTD  on,  Position  2) 
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A  separate  experiment  was  conducted  as  described  in  Knoll  to  determine  the  WE, 
W  D,  and  Wx  components  of  the  FWHM  (Knoll,  2000:  416-419).  The  FWHM  was 
determined  for  six  gamma  ray  energies  from  Co,  ~Eu  and  Cs  while  increasing  the 
detector  bias.  A  modification  to  one  CdTe  system,  as  shown  in  Figure  23,  was  performed 
to  facilitate  an  external  variable  power  supply  while  the  PX2T  performed  all  of  its  other 
power  and  amplification  tasks.  The  maximum  bias  voltages  applied  was  775  volts  based 
on  the  manufacturer’s  guidance  not  to  exceed  800  volts.  The  FWHM  was  also 
determined  with  an  ORTEC  precision  pulser.  FWHM  were  determined  for  pulser 
amplitudes  equivalent  to  100,  200  and  350  keV  for  each  bias  setting. 


Figure  23.  Modified  PX2T  for  external  power  supply. 

Next,  the  Nal(Tl)  detector  resolution  and  efficiency  were  evaluated  using  the 
geometry  shown  in  Figure  24.  Energy  spectra  were  collected  for  10  minutes  using  two 
sources,  137Cs  and  57Co.  A  sample  spectrum  is  shown  in  Figure  25.  The  resolution  was 
also  evaluated  at  four  positions  of  the  focus  adjustment  on  the  PA- 14  and  the  optimized 
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setting  was  recorded  and  used  throughout  the  remainder  of  the  thesis  research.  Efficiency 
calculations  included  peak  and  total  efficiencies  for  both  absolute  and  intrinsic 
measurements. 
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Figure  24.  Nal(Tl)  resolution  and  efficiency  geometry 


11200 
9600 
8000 
6400 
4800 
3200 
1600 
0 

0 


243 


486 

Energy(keV) 


729 


Figure  25.  Sample  137Cs  spectrum  using  Nal(Tl)  detector 
in  characterization  geometry 


System  Housing  Construction 

Alignment  of  the  source,  linear  collimator,  EEA,  and  conical  collimator  presented 
a  problem  to  establishing  the  Compton  spectrometer.  Independent  movement  of  the 
linear  collimator  and  the  EEA  relative  to  the  conical  collimator  is  needed  to  evaluate 
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alignment  and  to  provide  reproducible  geometries.  A  wooden  structure  shown  in  Figure 
26  was  created  to  facilitate  this  requirement.  The  structure  holds  the  PEA  flush  against  a 
Plexiglas  sheet  the  conical  collimator  rested  on.  This  maximizes  the  PEA  detector 
volume  in  line  with  the  cone  path.  The  linear  collimator  rests  on  a  Plexiglas  stage  over 
the  conical  collimator  and  EEA,  allowing  the  independent  adjustment  of  the  linear 
collimator  and  EEA  needed.  Component  alignment  techniques  included  the  use  of  a  laser 
and  a  solid  rod.  Both  were  directed  through  the  initial  collimator  for  alignment  with 
conical  collimator,  and  then  the  EEA  was  moved  into  position  without  disturbing  the 
collimators. 


Figure  26.  Compton  spectrometer  housing 
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Timing  Pulse  Alignment  and  Coincidence  Method  Evaluation 


Once  the  characterization  of  the  detectors  and  the  physical  geometry  of  the 
Compton  spectrometer  were  complete,  the  timing  electronics  were  established  and 
evaluated.  Both  the  T-SCA  and  PSA/T-SCA  were  used  to  produce  timing  pulses  from 
the  EEA  and  PEA  energy  pulses  and  tested  for  timing  variation  due  to  amplitude  walk.  A 
pulser  signal  was  sent  through  the  amplifiers  of  the  detectors  to  produce  simulated  energy 
signals  similar  to  those  created  when  a  gamma  ray  interacts  in  the  detectors.  Timing 
pulses  were  created  for  these  simulated  energy  pulses  using  the  T-SCA  and  PSA/T-SCA. 
Variations  in  the  timing  of  the  original  pulser  signal  and  the  timing  discriminator  signal 
were  measured  using  the  oscilloscope  measurement  tools  over  a  wide  range  of  pulser 
amplitudes.  Variations  were  attributed  to  the  effects  of  amplitude  walk. 

The  pulser  was  then  used  to  align  the  EEA  and  PEA  timing  pulses  with  each  other 
for  the  pulse  overlap  coincidence  method.  The  alignment  was  adjusted  using  the  built  in 
variable-delays  of  the  pulse  discriminators,  visually  established  using  the  oscilloscope, 
and  verified  by  counting  the  number  of  coincidence  events  over  a  set  time.  The  pulser 
operates  at  60  Hz  therefore  any  variation  in  that  rate  showed  the  coincidence  method 
operating  improperly.  The  PA- 14  did  not  provide  a  test  pulse  input,  and  the  pulser 
amplification  started  after  the  PEA’s  preamplifier.  A  Na  source  was  placed  between 
EEA  and  PEA  and  delay  curves  were  created  to  identify  the  optimal  pulse  alignment. 
Similarly,  the  pulser  was  used  to  test  the  coincidence  of  the  TAG  method.  The  output  of 
the  TAG  was  monitored  in  the  Genie  software  to  verify  proper  operation.  Again,  a  pulse 
rate  of  60  Hz  was  expected  under  the  prompt  coincidence  peak  and  variations  indicated 
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improper  operation.  Using  the  built  in  variable-delays  of  the  pulse  discriminators,  the 
shift  of  the  prompt  coincidence  peak  in  the  Genie  software  was  verified.  A  22Na  was  then 
placed  between  the  EEA  and  PEA  to  further  evaluate  the  TAG  method. 

Both  coincidence  methods  produce  a  logic  pulse  when  the  system  determines 
there  is  a  true  coincident  event.  The  triggering  of  this  logic  pulse  depends  on  the  method 
used  but  is  always  delayed  with  reference  to  the  energy  pulses.  The  energy  pulses  are 
sent  through  delay  amplifiers  for  alignment  with  the  logic  pulses  at  the  input  to  the 
ADCs.  The  ADCs,  when  operating  in  the  coincidence  mode,  process  only  the  energy 
pulses  that  have  a  corresponding  logic  pulse. 

Compton  Spectrometer  Data  Collection 

The  energy  spectra  from  the  EEA  and  PEA  were  collected  for  varying  conical 
collimator  gap  widths  using  both  coincidence  methods  and  both  initial  collimators. 
Because  the  spectra  were  expected  to  have  a  significant  dependence  on  the  alignment  of 
the  system  components,  graduated  scales  were  printed  and  affixed  to  the  EEA  and  the  top 
of  the  conical  collimator  as  shown  in  Appendix  E  and  Appendix  J.  The  graduated  scales 
were  used  to  determine  the  location  of  the  CdTe  crystal  in  the  detector  housing,  to 
evaluate  the  effects  of  the  RTD  on  the  active  detector  volume  and  to  reproduce 
component  alignments. 

The  location  of  the  CdTe  crystal  and  the  effects  of  the  RTD  were  detennined  by 

1  T7 

systematically  moving  the  EEA  through  a  collimated  Cs  beam  and  comparing  the 
count  rate  from  the  ICR  to  the  count  rate  from  the  RTD.  One -hour  counts  were 
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conducted  for  eight  locations  along  the  graduated  scale  on  the  EEA.  The  location 
corresponding  to  the  maximum  RTD  count  rate  was  used  as  the  alignment  point  for  the 
EEA  throughout  the  remainder  of  the  experiments. 

The  Compton  spectrometer  energy  spectra  from  the  EEA  and  PEA  were  collected 
for  ten  hours  using  137Cs.  The  alignment  of  the  initial  collimator  and  the  conical 
collimator  were  adjusted  and  analyzed  for  five  positions  using  the  graduated  scale  on  the 
conical  collimator.  The  alignment  process  was  done  in  two  steps.  First,  the  initial 
collimator  was  aligned  with  a  specific  grid  location  on  the  conical  collimator  graduated 
scale.  Then  the  EEA  was  moved  into  position  and  aligned  using  the  graduated  scale 
location  corresponding  to  the  maximum  RTD  count  rate.  Finally,  a  74-hour  count  was 
collected  using  22Na. 

Compton  Spectrometer  Simulation  Code 

A  program  was  written  to  simulate  the  Compton  spectrometer  in  order  to  assess 
the  effects  of  varying  geometric  parameters  on  the  resolution  and  efficiency  of  the 
system.  The  geometric  parameters  that  can  be  changed  include  the  dimensions  and 
relative  locations  of  the  initial  collimator,  the  EEA  detector  and  the  second  collimator. 
The  program  divides  the  source  into  equal  incremental  areas,  the  EEA  into  incremental 
volumes  and  the  gap  at  the  bottom  of  the  conical  collimator  into  equal  incremental  areas. 
This  creates  a  finite  number  of  three-“point’ ’  combinations  consisting  of  a  source  area,  an 
EEA  volume  and  a  conical  collimator  area.  All  possible  combinations  are  analyzed  by 
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performing  four  main  calculations  as  outline  below.  The  Mathematica  code  and  a 
detailed  description  of  its  function  are  given  in  Appendix  M. 

The  first  main  calculation  detennines  if  the  three-point  combination  is 
“geometrically  possible.”  The  program  uses  basic  geometric  and  algebraic  equations  to 
answer  two  questions.  First,  does  the  source  gamma  ray  pass  through  the  hole  at  the 
bottom  of  the  initial  collimator  before  reaching  the  EEA?  Second,  does  the  scattered 
gamma  from  the  EEA  pass  through  the  gap  of  the  conical  collimator  to  the  PEA?  If  both 
answers  are  positive,  the  three-point  combination  is  “geometrically  possible”  and  the 
remaining  four  calculations  are  performed. 

Since  the  source  gamma  passed  through  the  initial  collimator  hole  and  reached  the 
EEA,  the  next  question  is,  at  what  probability  does  the  isotropic  point  source  emit  gamma 
rays  into  the  top  surface  of  the  incremental  detector  volume?  This  probability  is 
described  as  the  solid  angle  subtended  by  the  incremental  detector  volume  from  the 
source  point.  The  program’s  second  calculation  uses  Equation  14  to  determine  this 
probability. 

The  source  gamma  has  now  reached  the  incremental  detector  volume  and  must 
Compton  scatter  at  a  specific  range  on  angles  to  pass  through  the  conical  collimator  gap. 
The  third  calculation  determines  the  probability  that  this  scatter  will  occur  in  two  steps. 
First  it  perfonns  the  solid  angle  integration  of  the  Klein-Nishina  fonnula,  Equation  4, 
over  the  specific  scatter  angles  needed  to  detennine  the  Compton  scatter  cross  section  for 
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that  three-point  combination.  This  cross  section  is  used  in  Equation  6  to  detennine  the 
probability  of  interaction  through  the  depth  of  the  incremental  detector  volume. 

The  fourth  calculation  “bins”  the  product  of  the  above  calculated  probabilities 
according  to  the  Compton  electron  energy  associated  with  the  three-point  combination. 
The  energy  is  found  by  using  the  average  scatter  angle  in  the  Compton  relation  described 
in  Equation  2.  The  program  continues  to  sum  the  probabilities  in  the  energy  channels 
until  all  three-point  combinations  have  been  examined.  The  result  is  an  energy  spectrum 
showing  the  probability  of  that  energy  occurring.  This  spectrum  can  then  be  used  to 
evaluate  different  designs  or  for  comparison  with  experimental  data. 
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V.  Data  Analysis  and  Results 


Chapter  Overview 

This  chapter  covers  the  experimental  and  simulated  data  collected  during  this 
research  and  is  divided  into  five  sections  following  the  methodology  of  Chapter  IV. 

1.  CdTe  and  Nal(Tl)  detector  characterization  including  resolution  and  efficiency. 

2.  The  effects  of  the  Compton  system  geometry  on  detector  efficiencies. 

3.  Analysis  of  the  start-stop  and  pulse  overlap  coincidence  timing  methods. 

4.  Analysis  of  experimental  energy  spectra  from  the  Compton  system. 

5.  Analysis  of  simulated  energy  spectra. 

CdTe  Characterization 

The  experimental  procedures  began  with  the  evaluation  of  the  CdTe  detector 
resolution  and  efficiency.  Eight  energy  spectra  were  taken  using  the  geometry  described 
in  Figure  20  and  are  located  in  Appendix  N.  The  resolution  data  is  summarized  in  Table 
2  and  Figure  27.  Generally,  the  detector  performed  above  the  expected  standards  given 
by  the  manufacturer  and  I  am  confident  it  is  operating  properly. 

Both  the  1.1  keV  FWHM  at  122  keV  and  the  5.27  keV  FWHM  at  662  keV  with 
the  RTD  active  were  better  than  the  Amptek’s  specification  of  1.239  keV  (Amptek,  2004: 
2)  and  5.9  keV  (Amptek,  2004:  28)  respectively.  No  specifications  were  given  for  the 
RTD  inactive,  but  Figure  18  shows  that  an  increase  in  FWHM  is  expected.  The  10.21 
keV  FWHM  at  662  keV  in  Position  2  with  the  RTD  active  is  larger  than  expected. 

Poorer  resolution  was  found  for  all  readings  from  Position  2  except  for  the  662  keV, 
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inactive  RTD  reading.  One  factor  for  this  general  trend  is  the  additional  scattering 
caused  by  the  nickel  detector  housing  in  Position  2  versus  the  beryllium  window  in 
Position  1 .  The  additional  scatters  create  a  distribution  of  photon  energies  which  causes 
an  increase  in  the  FWHM. 


Table  2.  CdTe  FWHM  Analysis  Data 


RTD 

Position 

Energy 

(keV) 

FWHM 

(keV) 

o(FWHM) 

(keV) 

Resolution 

(%) 

±o(Resolution) 

(%) 

On 

1 

122.0 

1.07 

0.01 

0.88 

0.01 

136.3 

0.98 

0.04 

0.72 

0.03 

661.7 

5.27 

0.15 

0.80 

0.02 

2 

121.8 

1.16 

0.02 

0.95 

0.02 

136.0 

1.16 

0.06 

0.85 

0.04 

662.7 

10.21 

0.22 

1.54 

0.03 

Off 

1 

121.9 

1.41 

0.03 

1.16 

0.02 

136.6 

1.13 

0.07 

0.83 

0.05 

661.4 

25.52 

1.5 

3.86 

0.23 

2 

122.4 

1.55 

0.04 

1.27 

0.03 

136.7 

1.4 

0.08 

1.02 

0.06 

661.6 

20.91 

0.5 

3.16 

0.08 

Quoted  resolution  from  Amptek 

RTD 

Position 

Energy 

(keV) 

FWHM 

(keV) 

Resolution 

(%) 

On 

1 

122.0 

1.24 

1.02 

On 

1 

662.0 

5.9 

0.89 
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Figure  27.  CdTe  resolution  vs.  gamma  ray  energy 

Using  the  same  energy  spectra  from  above  the  absolute  and  intrinsic  peak 
efficiencies  were  calculated.  The  absolute  peak  efficiencies  are  summarized  in  Table  3 
and  Figure  28  and  the  intrinsic  peak  efficiencies  are  summarized  in  Table  4  and  Figure 
29.  The  effects  of  the  RTD  on  and  off  and  in  Position  1  and  2  depend  on  the  energy  of 
the  gamma  ray.  As  expected  the  absolute  efficiency  for  the  low-energy  gamma  rays  was 
greater  in  Position  1  than  Position  2  with  the  RTD  both  on  and  off.  The  greater 
interaction  probability  associated  with  the  depth  advantage  achieved  in  Position  2  is  not 
significant  for  low  gamma  ray  energies.  This  is  further  demonstrated  by  the  fact  that 
there  is  only  a  fractional  increase  in  the  efficiency  with  the  RTD  switched  from  on  to  off 
for  low-energy  gamma  rays  entering  from  Position  1 .  The  high-energy  gamma  ray 
efficiencies  behave  differently.  High-energy  gamma  rays  entering  from  Position  1 
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experience  over  530%  increased  efficiency  with  the  RTD  switched  from  on  to  off.  The 
efficiency  with  the  RTD  on  also  increased  over  210%  from  Position  1  to  Position  2, 
demonstrating  the  increased  significance  of  the  depth  advantage  in  Position  2  with  higher 
energies. 


Table  3.  CdTe  Absolute  Peak  Efficiency  Data 


RTD 

Position 

Energy 

(keV) 

Net 

Peak 

(Counts) 

Net  Area 
Uncertainty 
(Counts) 

Absolute 

Efficiency 

(%) 

Absolute 

Eff. 

Uncertainty 

(%) 

On 

1 

122 

7740 

504.8 

2.08E-02 

6.68E-04 

136 

680 

57.8 

1.56E-02 

7.52E-04 

662 

15656 

153.4 

9.01E-05 

8.83E-07 

2 

122 

6720 

646 

9.23E-03 

2.99E-04 

136 

567 

73.8 

7.41E-03 

3.83E-04 

662 

9670 

105.4 

1.94E-04 

2.12E-06 

Off 

1 

122 

8170 

1056.9 

2.48E-02 

8.34E-04 

136 

744 

64 

1.66E-02 

1.04E-03 

662 

44729 

286.3 

4.80E-04 

3.07E-06 

2 

122 

7830 

917.3 

1.29E-02 

4.31E-04 

136 

849 

60.9 

9.96E-03 

5.71E-04 

662 

59802 

293 

3.23E-04 

1.58E-06 
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Figure  28.  CdTe  absolute  peak  efficiency  plot,  source-detector 
geometry  shown  in  Figure  20 

The  intrinsic  efficiency  calculations  take  into  consideration  the  solid  angle 
subtended  by  the  detector  and  the  decrease  in  detector  width  experienced  with  the  RTD 
active.  Position  2  intrinsic  efficiencies  are  greater  because  of  the  significant  decrease  in 
the  solid  angle  subtended  from  Position  1  (5mm  x  5mm)  to  Position  2  (1mm  x  5mm)  and 
the  increased  interaction  depth  of  Position  2.  An  estimated  decrease  in  detector  width  of 
70%  was  used  based  on  the  ratio  of  counts  from  the  RTD  to  the  ICR  found  when 
determining  the  location  of  the  crystal  inside  the  detector  housing.  The  highest  intrinsic 
efficiency  of  77%  was  found  from  Position  2  with  the  RTD  active. 
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Efficiency  (%) 


Table  4.  CdTe  Peak  Intrinsic  Efficiency  Data 


RTD 

Position 

Energy 

(keV) 

Solid  Angle 
(Steradians) 

Intrinsic 

Efficiency 

(%) 

Int.  Eff. 
Uncertainty 
(%) 

On 

1 

122 

0.0206 

12.688 

0.407 

136 

0.0206 

9.516 

0.459 

662 

0.0206 

0.055 

0.001 

2 

122 

0.0015 

77.325 

2.505 

136 

0.0015 

62.078 

3.209 

662 

0.0015 

1.625 

0.018 

Off 

1 

122 

0.0206 

15.128 

0.509 

136 

0.0206 

10.126 

0.634 

662 

0.0206 

0.293 

0.002 

2 

122 

0.0040 

40.527 

1.354 

136 

0.0040 

31.290 

1.794 

662 

0.0040 

1.015 

0.005 

Figure  29.  CdTe  peak  intrinsic  efficiency  plot,  source-detector 
geometry  shown  in  Figure  20 
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To  further  characterize  the  CdTe  detector  an  experiment  to  determine  the 
contributions  to  the  FWHM  was  conducted.  As  shown  in  Figure  30,  the  plot  of  the 

Pt-Wl  versus  the  inverse  high  voltage  showed  that  the  CdTe  does  not  behave  in  the 

same  way  the  Si(Li)  detector  described  by  Knoll  does  (Knoll,  2000:  418-419).  The 
expected  plot  would  show  a  decreasing  Wx  contribution  to  the  FWHM  as  the  high 
voltage  increased.  The  data  then  could  be  extrapolated  to  an  infinite  voltage  where  the 
contributions  of  Wx  would  be  minimal.  The  remaining  contribution  would  be 
contributed  only  to  WD.  The  y-intercept  of  the  extrapolated  lines  representing  the 
expected  contributions  of  WD  to  the  FWHM  and  an  estimated  Fano  factor  are  shown  in 
Table  5.  A  large  range  of  Fano  factors  from  0.007  to  1.61  was  found.  Fano  factors  range 
from  greater  than  zero  to  one  and  are  typically  less  than  1/3  for  semiconductors.  Only  the 
Fano  factor  for  the  59  keV  gamma  ray  of  0.25  was  of  the  expected  order  of  magnitude. 
Because  the  HVPS  was  not  increased  past  the  775  volts,  based  on  the  manufacturer’s 
warning,  it  is  unclear  if  a  high  enough  voltage  was  not  reached  or  if  as  Knoll  warns,  the 
carrier  velocity  became  saturated  at  the  higher  voltages  used.  Because  the  crystal 
dimensions  are  small  the  resulting  electric  field  is  approximately  7.75x1 05  volts/meter. 
Another  factor  that  affects  the  validity  of  this  method  is  the  pulse  selection  from  the 
RTD.  The  Fano  factor  describes  the  improvements  seen  from  expected  Poisson  statistics. 
If  the  RTD  is  preferentially  selecting  pulses  the  data  no  longer  represents  the  complete 
system  and  the  statistics  are  affected.  The  RTD  selection  is  expected  to  give  an  even 
greater  improvement  from  the  Poisson  statistics  resulting  in  a  lower  Fano  factor. 
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Figure  30.  Inverse  high  voltage  bias  versus  FWHM  plot 
Table  5.  Calculated  Fano  Factors  from  Figure  30 


Energy 

(keV) 

Y-Intercept 

Fano  Factor 

59 

0.60 

0.25 

88 

-1.11 

NA 

122 

1.48 

0.73 

244 

1.42 

0.34 

344 

0.24 

0.007 

662 

5.12 

1.61 

56 


Knoll’s  description  of  the  contributions  of  the  FWHM  for  a  lithium-drifted  silicon 
detector,  Si(Li),  (Knoll,  2000:  466-467)  was  also  used  to  describe  the  CdTe  detector 
FWHM  contributions.  A  linear  response  of  the  WT  -  We  for  gamma  energies  from  59 
keV  to  344  keV  suggested  that  the  contribution  of  WD“  was  dominant  in  this  range.  Fano 
factors  were  calculated  for  each  series  resulting  in  an  average  Fano  factor  of  0.205.  The 
resulting  contribution  of  Wx~=  WT“-  We  -Wd  is  shown  in  Figure  31.  The  leakage 
current  and  incomplete  charge  collection  contribution  to  the  FWHM  at  662  keV  is 
approximately  2.8  keV. 


CdTe  FWHM  Contributions 


Figure  31.  Estimate  of  the  FWHM  contributions  for  the  Amptek  CdTe  detector 
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Nal(Tl)  Characterization 


The  Nal  detector  performance  exceeded  the  manufacturer’s  specification.  The 
resolution  data  is  shown  in  Table  6  and  Figure  32.  A  resolution  of  6.7%  at  662  keV  was 
achieved  using  multinuclide  source  T108.  The  energy  spectrum  is  located  in  Appendix 
O. 


Table  6.  Nal(Tl)  Resolution  Data 


Nuclide 

Energy 

FWHM 

o(FWHM) 

R 

(keV) 

(keV) 

(keV) 

(%) 

Co57 

122.80 

11.01 

0.23 

9.0 

Cel39 

165.91 

15.08 

0.44 

9.1 

Hg203 

282.85 

24.00 

1.01 

8.5 

Snl  13 

382.18 

30.19 

1.51 

7.9 

Sr85 

514.47 

37.29 

0.71 

7.2 

Csl37 

661.76 

44.18 

0.18 

6.7 
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Nal(TI)  Detector  Resolution 


Figure  32.  Nal(Tl)  FWHM  versus  gamma  ray  energy 

For  an  accurate  measurement  of  the  absolute  and  intrinsic  efficiencies,  the  nearly 
mono-energetic  sources  of  57Co  and  I37Cs  were  used  in  the  geometry  shown  in  Figure  24. 
The  absolute  efficiency  data  is  found  in  Table  7.  The  intrinsic  efficiency  data  is  found  in 
Table  8.  The  intrinsic  total  efficiencies  of  96.5%  and  59.6%  for  122  keV  and  662  keV 
gamma  rays  respectively,  were  lower  than  shown  in  Knoll  of  100%  and  86%  (Knoll, 
2000:  337).  These  values  are  very  dependent  upon  source-detector  geometry  and  it  is 
suspected  that  the  values  from  Knoll  consider  additional  factors  such  as  the  detector 
window  material  and  source  self-attenuation  increasing  the  reported  efficiencies.  The 
662  keV  peak-to-total  ratio  of  0.42  was  found  to  be  reasonable  compared  to  a  reference 
value  of  0.5  (Knoll,  2000:  338),  further  confirming  the  proper  operation  of  the  detector. 
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Table  7.  Nal(Tl)  Absolute  Efficiency  Data 


Energy 

(keV) 

Total  Absolute 
Efficiency 
(%) 

Total  Abs.  Eff. 
Uncertainty 
(%) 

Total  Peak 
Efficiency 
(%) 

Total  Peak  Eff. 
Uncertainty 
(%) 

122 

1.00 

3.2E-05 

0.49 

1.5E-05 

662 

0.62 

1.5E-06 

0.26 

6.2E-07 

Table  8.  Nal(Tl)  Intrinsic  Efficiency  Data 


Energy 

(keV) 

Total  Intrinsic 
Efficiency 
(%) 

Total  Int.  Eff. 
Uncertainty 
(%) 

Peak  Intrinsic 
Efficiency 
(%) 

Peak  Int.  Eff. 
Uncertainty 
(%) 

122 

96.46 

3.1E-03 

46.89 

1.5E-03 

662 

59.56 

1.4E-04 

24.99 

4.1E-05 

Compton  System  Geometry  Effects 


The  system  was  setup  in  the  same  configuration  as  used  by  Capt.  Chris  Williams 
(Williams,  2003:  117).  A  complete  system  drawing  and  labeled  photograph  is  found  in 
Appendix  P.  This  section  analyzes  the  effects  of  the  geometry  and  clarifies  any  changes 
from  the  previous  experimental  work. 


The  effects  of  the  initial  collimator  on  the  system  resolution  have  been  found  to 
be  significant  and  required  detailed  analysis.  As  shown  in  Figure  33,  the  1.5mm  radius 
collimator  solid  angle  opens  up  at  the  detector  to  a  field  of  view  larger  than  the  entire 
detector.  It  is  important  to  understand  the  distribution  of  the  source  on  the  detector  based 
on  this  geometry.  The  initial  version  of  the  Mathematica  code  discussed  earlier,  created 
the  distribution  shown  in  Figure  34.  It  was  created  by  summing  the  solid  angle  subtended 
for  each  detector  point  by  all  points  on  the  source.  If  the  source  was  blocked  by  the 
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initial  collimator,  the  detector  point  was  not  seen  and  no  solid  angle  was  added.  As 
expected,  the  center  of  the  detector  is  seen  by  all  source  points  and  has  the  largest 
distribution  (white).  The  far  corners  of  the  detector  are  not  seen  by  the  entire  source  and 
therefore  have  the  smallest  distribution  (black).  This  shows  that  the  initial  collimator, 
although  not  ideal,  still  provides  some  reduction  in  the  angle  variance  between  the  source 
and  the  detector.  It  also  reduces  the  scattering  off  nearby  materials  into  the  detector 
which  add  to  the  background  noise.  A  smaller  collimator  hole  would  improve  the  desired 
alignment  of  the  source  and  detector  at  the  cost  of  further  reducing  efficiency. 


Figure  33.  Diagram  of  initial  circular  collimator  and  solid  angle 
subtended  by  CdTe  detector 
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Figure  34.  Source  distribution  on  CdTe  1  mm  x  5  mm  face  from  Position  2;  white 
pixels  represent  largest  distribution;  black  pixels  represent  smallest  distribution 

Reduction  in  the  conical  collimator  gap  also  affects  the  efficiency  of  the  system 
by  reducing  the  active  detector  volume  in  the  CdTe.  For  the  Compton  spectrometer 
system  to  work,  the  conical  collimator  must  be  able  to  select  only  those  gamma  rays 
scattered  at  the  desired  angle.  Using  a  similar  iterative  process  as  with  the  initial 
collimator  calculations,  the  program  analyzed  every  point  in  the  detector  to  determine  if  it 
passed  through  the  conical  collimator  gap  to  a  point  on  the  bottom  plane  of  the  collimator 
within  the  gap.  Figure  35  shows  the  three  dimensional  output  of  the  program 
representing  the  5x5x1  mm  detector  volume.  The  larger  points  represent  locations 
where  there  are  a  greater  number  of  possible  paths  through  collimator.  As  would  be 
expected,  a  region  in  the  detector  that  intersects  with  the  apex  of  the  collimator  gap 
would  have  the  most  possible  paths  through. 
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Figure  35.  Active  volume  of  CdTe  detector  based  on  conical  collimator  geometry; 
larger  circles  represent  a  larger  number  of  paths  through  the  conical  collimator  gap 
with  the  largest  occurring  at  the  apex  of  the  conical  collimator 

The  ratio  of  the  number  of  points  in  the  detector  volume  that  can  make  it  through 
the  conical  collimator  to  the  total  number  of  points  in  the  detector  is  a  good 
approximation  to  the  active  volume.  In  Figure  35  there  are  2336  points  that  make  it 
through,  out  of  the  3 125  total  points  giving  a  74.8%  usage  of  the  active  volume.  The 
more  refined  the  pixels  of  the  detector  and  the  pixels  of  the  collimator,  the  more  accurate 
the  calculation  will  be. 

Figure  35  suggests  that  the  alignment  of  the  conical  collimator  apex  with  the  EEA 
detector  volume  will  have  significant  affects  on  the  Compton  system.  In  order  to 
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improve  alignment  and  the  consistency  of  experimental  set  up,  the  location  of  the  CdTe 
crystal  inside  the  detector  housing  was  detennined  experimentally.  A  graduated  scale 
attached  to  the  detector  housing  was  used  to  produce  the  count  versus  position  graph 
shown  in  Figure  36.  The  detector,  with  the  RTD  on,  was  tested  at  eight  locations  each 
1/32  inch  (0.8  mm)  apart.  A  collimated  l37Cs  beam  with  a  3  mm  diameter  was  used. 
Based  on  the  data,  position  6  was  taken  as  the  detector  location  and  used  for  the 
remainder  of  the  experiment.  This  location  corresponds  to  1/8  inch  in  from  the  detector 
housing  face  and  is  in  agreement  with  Amptek’s  detector  dimensions  (Amptek,  2004: 
17). 


Figure  36.  CdTe  crystal  location  and  RTD  evaluation;  count  rate  from  RTD 
and  ICR  as  a  function  of  the  source  alignment  with  CdTe  graduated  scale 
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The  ratio  of  the  count  rate  from  the  RTD  and  the  ICR  was  also  used  to  evaluate 


the  effects  of  the  RTD  on  the  active  volume.  The  results  in  Figure  37  show  the  expected 
increase  in  the  RTD/ICR  ratio  as  the  source  is  moved  from  position  1  to  8.  As  the  source 
moves  to  the  front  edge  of  the  detector,  a  higher  count  rate  from  the  RTD  is  seen  as  the 
ICR  begins  to  decrease.  When  the  source  is  over  position  6  a  70%  reduction  in  detector 
volume  is  observed  with  the  RTD  active. 


Figure  37.  RTD/ICR  Ratio;  determination  of  active  RTD  on  detector  volume 
Coincidence  Electronics  Evaluation 

The  overlap  method  of  coincidence  was  setup  as  follows.  First,  SCAs  were  used 
to  produce  logic  pulses  from  the  CdTe  and  the  Nal(Tl)  energy  pulses  based  on  their 
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arrival  times  as  shown  Figure  38.  These  pulses  were  initially  aligned  visually  using  the 
oscilloscope  and  the  adjustable  delay  in  the  SCAs. 
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Figure  38.  Oscilloscope  output-SCA  logic  pulse  generated 
from  EEA  and  PEA  energy  pulses 

These  pulses  were  sent  to  GDGs  for  pulse  width,  height  and  delay  adjustment.  The 
GDGs  outputs  were  then  sent  to  the  LGSC.  The  LGSC  uses  the  pulse  overlap  method  as 
shown  in  Figure  39.  A  logic  pulse  is  produced  in  the  region  that  the  pulses  overlap. 
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Figure  39.  Oscilloscope  output-LGSC  logic  pulse  created  by  coincident 

EEA  and  PEA  timing  pulses 

The  LGSC  output  is  then  sent  to  the  gate  input  of  the  ADC.  The  ADC  produces  a  logic 
acceptance  pulse  when  it  receives  a  pulse  to  the  ADC  In.  In  the  delayed  mode,  the 
acceptance  gate  can  be  extended  past  the  end  of  the  pulse  going  to  the  ADC  In  as  shown 
in  Figure  40.  In  coincidence  mode  the  energy  pulse  is  accepted  when  the  gate  in  pulse 
overlaps  the  ADC  delayed  linear  gate. 
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Figure  40.  Oscilloscope  output-ADC  delayed  linear  gate  alignment  with  LGSC 

output 

The  visual  alignment  was  then  refined  by  creating  the  delay  curve  shown  in  Figure  41. 
The  delay  that  corresponds  to  the  highest  counts  is  the  optimized  solution.  The  very  low 
count  rate  seen  makes  the  optimization  less  accurate.  A  similar  process  is  used  to 
establish  the  start-stop  method. 
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Figure  41.  Delay  curve-counts  versus  adjustable  delay 
to  determine  optimized  alignment 

The  TAC  method  produces  a  pulse  with  a  voltage  proportional  to  the  timing 
between  the  EEA  and  PEA  pulses  that  can  be  gated  using  a  SCA  as  shown  in  Figure  42. 
The  SCA  output  represents  a  coincidence  event  and  is  sent  to  the  ADC  gate  in  as 
described  above.  Because  the  only  pulse  timing  requirement  is  the  arrival  of  the  stop 
pulse  (PEA)  after  the  start  pulse  (EEA),  the  TAC  method  does  not  require  exact 
alignment  of  pulses.  Both  overlap  and  start-stop  were  used  in  the  Compton  spectrometer 
with  similar  results. 
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Figure  42.  Oscilloscope  output  using  TAC  method;  time  between  CdTe  and  Nal(Tl) 

determines  TAC  output  pulse  height 

Experimental  Compton  Spectrometer  Energy  Spectra 

The  experimental  Compton  spectrometer  data  collected  was  not  the  same  as  found 
during  the  previous  research  of  Capt.  Williams.  Included  for  analysis  are  four  spectra;  a 

Cs  pulse-overlap  method,  a  Cs  start-stop  method,  a  Na  start-stop  method  and  the 

22 

previous  "Na  spectrum  “smoothed”,  shown  in  Figure  43,  Figure  44,  Figure  45  and  Figure 

i  -jn 

46,  respectively.  The  Cs  spectra  show  the  majority  of  counts  roughly  centered  about 
100  keV.  This  suggests  that,  to  some  extent,  the  second  collimator  is  preferential  to  the 
scatter  gamma  rays  corresponding  to  the  collimators  scatter  angle.  This  is  reinforced 
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with  the  "Na  spectrum  centered  roughly  around  61  keV,  the  expected  Compton  electron 
energy  from  the  551  keV  gamma  ray.  Although  the  spectra  suggest  that  the  counts 
collected  are  real,  these  spectra  are  in  stark  disagreement  with  the  previous  research 
which  showed  a  high  resolution  peak  near  the  expected  energy  of  100  keV.  Although  the 
cause  of  this  discrepancy  remains  unknown  computer  simulation  has  provided  some 
additional  information. 
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Figure  43.  10-Hour  Compton  spectrometer  spectrum-137Cs  source  using  pulse- 

overlap  method 
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Figure  44.  10-Hour  Compton  spectrometer  spectrum- 137Cs  source  using  start-stop 

method 


71 


40  80 

Energy(keV) 


Lett  Marker:  116  :  39.1  keV  FWHM.  FWTM:  8.983. 25.337  keV 

Right  Marker:  240  :  80.9  keV  Gaussian  Ratio:  1.547 


Centroid:  171  :  57.5keV  ROI  Type: 

Area:  1951  ±4.43%  Integral:  3076 


Figure  45.  74-Hour  Compton  spectrometer  spectrum-22Na  source 
using  start-stop  method 
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Figure  46.  13  point  smoothed  74-Hour  Compton  spectrometer 
spectrum-22Na  source  using  start-stop  method 


72 


Simulated  Compton  Spectrometer  Energy  Spectra 


The  Mathematica  code  written  to  simulate  the  operation  of  the  Compton 
spectrometer  was  executed  for  a  wide  variety  of  geometric  scenarios  for  verification  of 
the  code’s  perfonnance,  for  comparison  with  experimental  data  and  for  detennining  what 
geometries  would  meet  the  requirements  of  this  project.  The  geometric  parameters 
analyzed  for  effects  on  system  resolution  and  efficiency  were  the  initial  collimator  radius, 
the  conical  collimator  gap,  and  the  detector  dimensions  and  location  relative  to  the 
conical  collimator  apex. 

The  simulated  spectra  below  represent  counts  versus  gamma  ray  energy  for  1- 
hour  collections  using  the  combined  137Cs  source.  It  is  also  assumed  that  the  Nal(Tl)  has 
a  25%  peak  intrinsic  efficiency  which  represent  the  effects  of  using  an  energy  gate  on  the 
PEA  since  coincident  events  are  counted  only  if  the  scattered  gamma  ray  is  fully 
absorbed  in  the  PEA.  Although  many  plots  will  show  very  low  and  often  fractional 
counts,  the  1-hour  collection  time  was  used  throughout  to  facilitate  efficiency 
comparisons  from  one  system  geometry  to  the  next.  Because  the  effects  on  resolution 
and  efficiency  are  dependent  on  the  ratios  of  the  geometric  parameters,  the  comparisons 
below  are  used  to  evaluate  specific  system  geometries  and  attempts  only  very  broad 
classification  of  these  effects. 

First,  the  simulation  verified  the  expected  reduction  in  the  FWHM  and  efficiency 
associated  with  reducing  the  conical  collimator  gap.  The  simulation  predicted  a  19% 
reduction  in  FWHM  with  an  81%  reduction  in  efficiency  with  a  reduction  in  the  gap  from 
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1.0  mm  to  0.5  mm  as  shown  in  Figure  47  and  Figure  48.  This  disproportionate  reduction 
shows  that  the  high  concentration  of  available  paths  through  the  conical  collimator  from 
the  apex  area  is  the  main  contribution  to  variance  in  angles  and  therefore  to  FWHM. 

This  conclusion  is  further  substantiated  when  evaluating  the  variation  in  detector  location 
relative  to  the  conical  collimator  apex.  Next,  the  simulated  reduction  of  the  initial 
collimator  radius  from  1.0  mm  to  0.25  mm  resulted  in  a  50%  reduction  in  FWHM  and 
61%  in  efficiency.  Shown  in  Figure  49,  this  significant  reduction  in  initial  radius  changes 
the  distribution  in  Figure  34  to  a  small  region  of  high  distribution  at  the  center  of  the 
detector  face.  The  predicted  FWHM  of  4.6  keV  demonstrates  that  reducing  the  initial 
collimator  alone  will  not  provide  the  resolution  needed.  The  surprising  result  of  an 
increased  FWHM  was  found  when  the  detector  width  was  reduced  from  0.5  mm  to  0.3 
mm.  The  expected  result  was  a  reduction  in  the  FWHM  because  of  a  decreased  scatter 
angle  distribution  and  a  reduced  efficiency  because  of  a  reduced  detector  volume.  The 
small  increase  in  the  FWHM  shown  in  the  comparison  of  Figure  47  and  Figure  50  is  a 
result  of  a  larger  relative  reduction  in  efficiency,  and  consequently  a  lower  half 
maximum,  without  a  significant  reduction  in  angle  distribution.  This  effect  demonstrates 
the  significant  dependence  on  the  ratios  of  the  geometric  parameters  as  this  effect  would 
not  be  seen  with  an  initial  collimator  of  a  smaller  radius. 
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Figure  47.  Simulated  Compton 
spectrometer  spectrum  -  conical 
collimator  gap=1.0  mm 


Figure  49.  Simulated  Compton 
spectrometer  spectrum  -  collimator 
radius=0.25  mm 
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Figure  48.  Simulated  Compton 
spectrometer  spectrum  -  conical 
collimator  gap=0.5  mm 
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Figure  50.  Simulated  Compton 
spectrometer  spectrum  -  detector 
width=0.3  mm 


The  effects  of  detector  location  relative  to  the  conical  collimator  apex  is  explored 
in  Figure  51,  Figure  52,  and  Figure  53;  representing  the  apex  just  above  the  top  of  the 
detector,  in  the  center  of  the  detector  and  just  below  detector,  respectively.  Moving  the 
apex  out  of  the  detector  volume  decreased  the  FWHM  by  50%  for  both  above  and  below 
the  detector.  As  mentioned  above,  the  high  concentration  of  available  paths  through  the 
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conical  collimator  from  the  apex  area  is  the  main  contribution  to  the  FWHM.  As  the 
center  of  the  “X”  distribution,  shown  in  Figure  35,  moves  above  or  below  the  detector, 
respectively  only  the  “legs”  or  “arms”  of  the  distribution  remain  in  the  detector  volume. 
These  portions  of  the  distribution  have  a  smaller  angle  variance  because  they  are  limited 
to  scattering  through  only  a  portion  of  the  collimator.  The  scatter  angles  associated  with 
the  “legs”  are  smaller  and  as  expected  the  peak  in  Figure  5 1  is  shifted  below  100  keV. 
The  “arms”,  on  the  other  hand,  have  larger  scatter  angles  and  the  peak  in  Figure  53  is 
shifted  above  100  keV.  The  predicted  results  with  the  apex  above  and  below  the  detector 
showed  similar  reductions  in  FWHM  and  efficiency.  Therefore  moving  the  apex  above 
the  detector  would  be  preferred  because  of  the  improved  resolution  associated  with 
collecting  the  lower-energy  Compton  electrons. 
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Figure  51.  Simulated  Compton  spectrometer  spectrum  -  detector  location=2.0  mm 
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52.  Simulated  Compton  spectrometer  spectrum  -  detector  location=-1.0  mm 


Figure  53.  Simulated  Compton  spectrometer  -  detector  location=-4.0  mm 
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Using  the  geometry  of  the  experimental  system,  a  FWHM  of  13.0  keV  was  found 
as  shown  in  Figure  54.  This  value  is  approximately  15  keV  better  than  an  estimated 
experimental  value  from  Figure  43  and  12  keV  worse  than  previous  research  values.  The 
simulated  one  hour  spectrum  using  the  combined  137Cs  source  predicted  130  counts  under 
the  curve.  Current  experimental  data  had  305  counts  under  the  curve  in  a  10-hour 
spectrum.  Previous  experimental  data  had  435  counts  under  the  curve  in  a  1-hour 
spectrum. 
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Figure  54.  Simulated  1-hour  Compton  system  spectrum  using  combined  137Cs 
source  showing  counts  versus  energy 

The  program  was  also  used  to  produce  a  1-hour  '  Na  spectrum  shown  in  Figure 
55  for  comparison  with  experimental  results  shown  in  Figure  45  and  Figure  46.  First  it  is 
important  to  note  that  the  simulation  demonstrated  the  natural  improvement  in  resolution 
due  to  the  reduced  gamma  ray  energy  of  5 1 1  keV.  For  example,  using  a  range  of  0  as 
0.5 1  to  0.55  radians  in  Equation  2,  the  Compton  electron  energy  has  a  range  of  12.5  keV 
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for  an  incident  gamma  ray  of  662  keV  and  only  8.0  keV  for  a  5 1 1  keV  incident  gamma 
ray.  The  Genie  software  found  a  9  keV  FWHM  for  the  unsmoothed  spectrum  and  13.6 
keV  for  the  13-point  smoothed  spectrum  compared  to  the  simulation  FWHM  of  8.3  keV. 
The  small  improvement  in  agreement  is  attributed  to  a  higher  activity  source,  longer 
collection  time  and  tighter  energy  requirements  in  the  PEA.  There  remains  a  strong 
disagreement  with  predicted  and  measured  efficiencies.  The  total  count  measured  in  a 
74-hour  period  was  3106,  while  the  simulation  predicted  23236  in  the  same  amount  of 
time.  The  disagreement  in  efficiency  is  because  the  simulation  does  not  take  into 
consideration  the  reduced  efficiency  of  the  Nal(Tl)  in  the  Compton  spectrometer 
configuration. 
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Figure  55.  Simulated  1-hour  Compton  system  spectrum  using  60  pcurie  22Na  source 

showing  counts  versus  energy 

In  the  Compton  spectrometer  configuration  only  a  small  portion  of  the  Nal(Tl) 
crystal  is  intersected  by  the  path  through  the  conical  collimator.  In  the  characterization  of 
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the  Nal(Tl)  detector  a  peak  efficiency  of  25%  was  found  using  the  entire  detector  volume 
with  the  incident  gamma  rays  nearly  normal  to  the  detector  face.  The  resulting  linear 
distance  through  the  detector  is  approximately  three  inches.  In  the  Compton 
spectrometer  configuration  the  gamma  rays  enter  the  detector  at  a  30°  angle  resulting  in  a 
one  inch  linear  distance  through  the  crystal.  Because  peak  events  correspond  not  only  to 
direct  photoelectric  absorption  events  but  also  Compton  events  followed  by  photoelectric 
absorption,  it  is  difficult  to  accurately  estimate  the  reduction  in  efficiency  based  on  the 
linear  distance  through  the  crystal.  The  reduction  of  Compton  events  followed  by 
photoelectric  absorption  will  be  greater  than  the  reduction  in  the  direct  photoelectric 
events;  therefore  a  one  third  reduction  is  conservative.  This  efficiency  correction  would 
decrease  the  simulation  prediction  of  23236  counts  to  less  than  7700  counts. 

Although  simulation  and  experimental  results  are  not  explicitly  linked,  there 
remains  value  in  using  the  program  to  detennine  what  geometric  parameters  are  needed 
to  achieve  the  required  resolution  set  forth  at  the  start  of  this  research.  A  FWHM  of 
0.533  keV  was  achieved  for  a  system  geometry  of  0.25  mm  initial  collimator  radius,  0.10 
mm  conical  collimator  gap,  detector  dimensions  of  4.5  x  0.3  x  4.5  mm  ,  and  raising  the 
conical  collimator  2.5  mm  closer  to  the  detector.  The  spectrum,  shown  in  Figure  56, 
required  a  200  curie  source  counted  for  one  hour  to  acquire  150  counts  under  the  curve. 
This  geometry  uses  the  “legs”  of  the  detector  volume  “X”  distribution  and  therefore  the 
peak  is  shifted  down  to  approximately  96  keV. 
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Figure  56.  Simulated  662  keV  spectrum  with  required  resolution 
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VI.  Conclusions  and  Recommendations 


Conclusions  of  Research 

Three  investigated  questions  were  used  to  focus  the  research  of  this  thesis  and  are 
answered  here. 

1 .  Can  the  previous  results  be  reproduced?  I  was  not  able  to  reproduce  the  results  of 
the  previous  research.  This  clearly  marks  these  finding  for  skepticism  and  raises  many 
more  questions.  The  important  question  is  which  results  more  closely  represent  the 
actual  characteristics  of  the  detector  system?  I  start  my  argument  in  support  of  my  results 
based  on  the  simple  fact  of  geometric  calculations.  The  rough,  two-dimensional 
calculations  of  the  system  produce  a  FWHM  of  approximately  14  keV.  This  calculation 
does  not  take  into  consideration  the  contributions  to  FWHM  from  the  detector  and 
electronics  and  therefore  is  a  best  case  scenario.  My  experimental  results,  although  not 
decisive,  have  general  characteristics  that  are  similar  to  those  expected.  The  previous 
results  are  significantly  better  than  expectations.  There  are  some  possible  solutions  for 
this  improvement  such  as  a  non-perpendicular  alignment  of  the  initial  collimator  to  the 
conical  collimator  or  an  improved  collimation  of  the  source  due  to  misalignments  of  the 
source,  initial  collimator  and/or  detector  volume.  Currently  no  reason  has  been  verified. 

2.  How  does  the  energy  gating  affect  the  Compton  Spectrometer?  The  main  effect 
of  adding  the  energy  gate  to  the  PEA  is  a  further  reduction  in  efficiency.  The  Nal(Tl) 
detector  is  chosen  for  its  very  high  efficiency  to  ensure  that  a  minimal  number  of  true 
coincident  events  are  lost  because  of  failure  of  the  scattered  gamma  ray  to  interact  in  the 
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PEA.  The  purpose  of  the  energy  gate  is  to  suppress  events  that  are  due  to  other  sources 
outside  the  energy  window  of  interest  and  to  further  reduce  accidental  coincidence.  This 
requires  a  better  knowledge  of  the  scattered  gamma  ray  energy.  The  energy  gate  selects  a 
window  around  the  range  of  full-energy  peaks  of  the  scattered  gamma  rays  of  interest. 
The  efficiency  of  the  Nal(Tl)  is  reduced  from  its  total  intrinsic  efficiency  to  the  peak 
intrinsic  efficiency.  For  662  keV  gamma  rays  this  was  found  to  be  a  reduction  from  60% 
to  25%.  The  efficiency  is  further  reduced  by  a  factor  of  at  least  1/3  because  of  the 
reduction  in  the  amount  of  Nal(Tl)  crystal  intersecting  the  path  of  the  conical  collimator. 

3.  Can  modeling  the  Compton  spectrometer  help  with  future  design?  Yes,  modeling 
of  the  system  is  useful  for  evaluating  new  system  designs  before  construction.  The 
interdependence  of  the  geometric  parameters  is  not  always  clear  and  modeling  can 
identify  unexpected  geometric  dependence  before  new  systems  are  constructed.  With 
specific  resolution  and  efficiency  requirements  a  system  design  can  be  evaluated  for 
feasibility. 

Recommendations  for  Action 

Research  of  the  Compton  spectrometer  should  continue  for  several  reasons.  First, 
real  world  applications  require  the  type  of  resolution  this  thesis  tried  to  achieve.  The 
concepts  of  the  system  are  well  understood  and  there  are  many  alternative  designs  that 
can  be  researched  to  overcome  the  major  disadvantage  of  efficiency.  Second,  the 
discrepancy  with  previous  experimental  data  must  be  further  researched  to  truly  evaluate 
the  Compton  spectrometer  system.  Lastly  is  the  educational  benefit  for  future  students. 
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This  thesis  covered  a  large  portion  of  the  nuclear  engineering  curriculum  and  was  a  great 
way  to  tie  this  educational  experience  together. 

Recommendations  for  Future  Research 

As  mentioned  above,  the  major  disadvantage  of  the  system  is  efficiency.  Further 
research  should  focus  on  overcoming  this  shortcoming.  There  are  several  possible 
avenues  that  can  be  evaluated.  The  previous  research  of  Capt.  Williams  used  a  two-angle 
collimator  to  increase  efficiency  while  maintaining  only  one  PEA  to  collect  the  scattered 
gamma  rays.  Multiplexing  the  scattered  angles  is  essential  to  increasing  efficiency,  but 
requires  a  method  to  distinguish  absorption  in  each  angle  separately.  Conical  collimators 
can  be  created  with  scintillation  material  with  each  conical  volume  separated  by  a 
reflective  material. 

The  research  should  also  include  improving  the  computer  simulation.  There  are 
two  main  areas  that  can  be  further  developed.  First,  the  fidelity  of  the  simulation  is 
currently  limited  to  the  selected  “pixel  size”  of  each  component.  This  limitation  can  be 
reduced  by  writing  the  program  in  FORTRAN  or  other  similar  program  with  a  more 
advanced  compiling  system.  The  simulation  could  also  be  easily  adapted  to  a  Monte- 
Carlo  method  that  would  not  have  the  dependence  on  “pixel  size.”  Second,  the  program 
requires  the  flexibility  to  evaluate  off-alignment  system  configurations.  Currently,  the 
program  is  written  so  that  changes  to  geometric  parameters  do  not  affect  the  vertical 
alignment  of  the  system.  Off-alignment  simulations  would  be  valuable  to  assess  the 
sensitivity  of  that  parameter  and  detennine  design  specifications. 
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Appendix  A:  239Pu  and  240Pu  Spectra 


The  following  energy  spectra  demonstrate  the  complex  characteristics  of 
plutonium  samples.  The  plots  were  produced  by  ht  Idaho  National  Engineering  and 
Environmental  Gamma  Ray  Spectrometry  Center’s  Online  Catalogue  (IINEE). 
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Figure  57.  239Pu  Energy  spectrum  0-400  keV,  (INEE,1999) 
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Figure  58.  239Pu  Energy  spectrum,  400-800  keV  (INEE,1999) 
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Figure  59.  240Pu  Energy  spectrum,  0-800  keV  (INEE,1999) 
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Appendix  B:  Summary  of  Previous  Experimental  Findings 


The  following  provides  a  summary  of  Capt.  Williams’  experimental  findings 
published  in  Plutonium  Isotopic  Ratio  Determination  Using  Compton  Spectrometer 
System. 
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Figure  60.  One-hour  single  collimator  Compton  spectrometer  spectrum,  channel 
width  0.050±0.002  inch  (Williams,  2003:  80) 


Table  9.  Spectral  data  from  Compton  spectrometer  (Williams,  2003:  81) 


Collimator  Geometry 

Max 

Counts 

Channel 

Number 

Energy 
in  keV 

Area 

Counts 

Integral 

Counts 

FWHM 
in  keV 

nVTM 

in  keV 

Single 

Angle 

Channel 

in  inches 

0.05<K002 

43 

6526 

100.753 

481 

23.7 
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0.194 

0.629 

0.050^.002 

48 

6530 
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24.5 
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44 
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100.734 
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- 

25.5 
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38 

6535 

100.928 
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- 

23.5 

491 
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0.638 

0. 040^002 

32 

6530 

100.831 
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21.1 

432 
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37 

6536 

:c:.94 
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22.1 
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0.644 
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22.3 
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33 
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- 
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31 
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- 

20.2 
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0.638 
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0.02CK002 

31 

6527 

100.773 

255 

23.5 

301 

0.174 
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Figure  61.  Plot  of  average  FWHM  from  100  keV  peak  using  Compton  spectrometer 

(Williams,  2003:  85) 


Figure  62.  Simulated  wgPu  spectrum  using  Compton  spectrometer 

(Williams,  2003:  96) 
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Figure  63.  Simulated  rgPu  spectrum  using  Compton  spectrometer 

(Williams,  2003:  96) 
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Appendix  C:  Klein-Nishina  Approximation  Evaluation 


The  Klein-Nishina  formula  is  used  in  the  system  simulation  code  to  determine  scatter 
probabilities  inside  the  EEA;  therefore  an  understanding  of  the  limitations  of  this 
approximation  is  necessary.  Comparison  calculations  of  the  Klein-Nishina  formula  to  an 
empirical  formula  (Massaro  and  Matt,  1986:  545-547)  and  to  XMuDat  (Nowotny,  1998) 
follow. 


Klein-Nishina  calculations: 


Acd=  112.411  ; 
ATe  =  127 . 60  ; 


ZcdTe  =  100  ; 

PcdTe  =6.06  (*Density  from  Knoll  *)  ; 

6.02214  xlO 23 
NcdTe  =  — 77 - 7 — ; —  ' 

(■Aa  +  ATe  ) 

ro  =  100  *  2.817940285  *  10~15  ; 

KN  re  ,  Z  ,  hv  1  = 

_ q[hv]2  (1-Cos  [e])2 _ 

,  (1  +Cos  [0]2)  (1  +  a[hv]  (1  -Cos  [6]  )) 


f  1  f 

'  1  +  Cos  [0] 2  ' 

l  1  +a[hv]  (1  -  Cos  [0])  ) 

2 

CdTe  KNDat  =  TableFozm  [Table  [{x,  NcdTe  2  7r  NIntegrate  [KN[0,  ZcdTe,  x]  *Sin  [6]  ,  {©,  0,  7r}]}, 
{x,  10,  1000  ,  10  }],  TableHeadings  {None  ,  {"y  (keV  )  "  ,  "  o(cm2/g)  "  }  }  ] 


Y  (keV) 

o (cm2 /g) 

100 

0.123636 

200 

0.101991 

300 

0.0886846 

400 

0.0794585 

500 

0.0725551 

600 

0.0671181 

700 

0.062678 
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Empirical  formula  calculations  (Massaro  and  Matt,  1986:  545-547): 


p[Z_]  =  Z  .6652  (*10-28*)  ; 

ai  =  2.0603  ; 

a2  =  1.1691  *  10  _1 ; 

bi  =  5.9973  ; 

b2  =  2.9267  ; 

b3  =  8.2038  *  10  “2; 

61.74  +  65.174 

cx  =  - - - ;  (*cd  :  61.714  ,  Te  :  65.174  *) 

7.7822  *102  +  5.5160  *  102  ,  , 

c2  =  - ;  (*  Cd:  7.7822  *10 2 ,  Te  :  5.5160  *102  *) 

2. 

32.736  +  35.66 

K  =  - ;  (*  Cd:  32.736  ,  Te :  35.660  *) 

2. 


1  +  ai  EE  +  a2  EE2  K 

cs  [EE_  ,  Z_]  =  p  [ Z]  - - - - -  -  - - -  ; 

1  +  bi  EE  +  b2  EE^  +  b3  EE  J  1  +  ci  EE  +  c2  EEZ 


Newi^jp  =  TableForm  [Table  [|l0  x,  2.  Nphto  10~24  cs  [  — ,  50.18  ]},  {x,  1,  70,  10}], 

TableHeadings  (None,  {  "y  (keV )  "  ,  "ct  (an2/g)  "  }  }  ] 
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Figure  64.  Comparison  plot  of  Klein-Nishina  formula,  an  empirical  formula 
(Massaro  and  Matt,  1986:  545-547)  and  XMuDat  (Nowotny,  1998) 
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Appendix  D:  Compton  Relation  Calculations 


The  following  calculations,  produced  in  Mathematica®  demonstrate  the  relation 
between  the  resolution  achieved  in  collecting  the  Compton  electron  in  the  EEA  and  the 
calculated  gamma  ray  full-energy  peak  resolution.  The  calculations  are  performed  for 
200-700  keV  gamma  rays  with  Compton  scatter  angles  of  0.4,  0.53,  0.6,  and  0.7  radians 


Formulas  used: 


mo  =  511; 

x 

CarptElectre  ,  x  1  ;  =  x - ; 

-  -  1+  x  (1  -  Cos[0] ) 

1T[0  ,  CE  ]  =  - 1 -  (CE-  CECos[0]  -VCE  V-f+C os[0]  V^CE4CECos[0]  -4  mo); 

-  -  2(1-  Cos[0] )  v  ' 


CcrnptEWHM  [  Newe_ ,  EWHM_]  = 

Re[y[New0,  |carptElect[New0,  r] 


“)]] 


Re[y[New0,  |corrptElect[New0,  r] 


FWHM 

2. 


Table  10.  Compton  electron  energy  as  a  function  of  incident  gamma  ray  energy  and 

scatter  gamma  ray  angle  (in  radians) 


Incident 

Gamma 

Ray 

(keV) 

Compton 

Electron 

(keV) 

0=4 

Compton 

Electron 

(keV) 

0=53 

Compton 

Electron 

(keV) 

0=6 

Compton 

Electron 

(keV) 

0=7 

200 

6.0 

10.2 

12.8 

16.9 

300 

13.3 

22.4 

27.9 

36.4 

400 

23.3 

38.8 

48.1 

62.2 

500 

35.9 

59.2 

73.0 

93.5 

600 

50.9 

83.2 

102.1 

129.8 

700 

68.3 

110.7 

135.1 

170.6 
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Table  11.  Compton  Electron  FWHM  vs.  Incident  Gamma  Ray  FWHM  (y=200  keV) 


EEA 

FEP 

FEP 

FEP 

FEP 

FWHM 

FWHM 

FWHM 

FWHM 

FWHM 

(keV) 

(0=.4) 

(0=  .53) 

(0-  •  6) 

(0-  •  7  ) 

0 . 1 

1  .694 

1  .007 

0 .8072 

0.6194 

0.2 

3.388 

2 .014 

1  .  614 

1.239 

0 . 3 

5 .082 

3.021 

2 .422 

1 .858 

0 . 4 

6.776 

4  .  027 

3.229 

2 .477 

0 . 5 

8.470 

5 .034 

4 .036 

3 .097 

0  .  6 

10.17 

6 .042 

4  .844 

3.716 

0 . 7 

11.86 

7.049 

5 . 651 

4 .336 

00 

o 

13.56 

8 .056 

6.459 

4  .  955 

0  .  9 

15.25 

9.064 

7.266 

5 .575 

1  . 

16.95 

10.07 

8  .074 

6.194 

1  . 1 

18.65 

11.08 

8  .881 

6.814 

1 . 2 

20.35 

12.09 

9.689 

7 .433 

1 . 3 

22 . 05 

13.10 

10.50 

8 . 053 

1  .  4 

23.75 

14.10 

11.31 

8  .  673 

Table  12.  Compton  Electron  FWHM  vs.  Incident  Gamma  Ray  FWHM  (y=300  keV) 


EEA 

FEP 

FEP 

FEP 

FEP 

FWHM 

FWHM 

FWHM 

FWHM 

FWHM 

(keV) 

(0- • 4 ) 

(0=  .53) 

(0=  •  6) 

(0=  •  7 ) 

0 . 1 

1 .154 

0 . 6967 

0.5638 

0 .4388 

0.2 

2.309 

1.393 

1.128 

0 .8776 

0 . 3 

3.463 

2 .090 

1  .691 

1.316 

0 . 4 

4.618 

2.787 

2 .255 

1 .755 

0 . 5 

5 .773 

3.484 

2.819 

2.194 

0  .  6 

6 . 927 

4.181 

3.383 

2 . 633 

0 . 7 

8 .082 

4  .877 

3 . 947 

3 . 072 

o 

00 

9.237 

5 .574 

4.511 

3.510 

0  .  9 

10.39 

6.271 

5 .075 

3.949 

1 . 

11.55 

6  .  968 

5 . 638 

4 .388 

1 . 1 

12.70 

7 .665 

6.202 

4 . 827 

1 . 2 

13.86 

8.362 

6.766 

5.266 

1 . 3 

15.01 

9.059 

7.330 

5.704 

1  .  4 

16.17 

9.756 

7.894 

6 .143 
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Table  13.  Compton  Electron  FWHM  vs.  Incident  Gamma  Ray  FWHM  (y=400  keV) 


EEA 

FEP 

FEP 

FEP 

FEP 

FWHM 

FWHM 

FWHM 

FWHM 

FWHM 

(keV) 

(0=.4) 

(0= ,  5  3) 

(0-  •  6) 

(0-  •  7  ) 

0 . 1 

0.8849 

0.5419 

0 .4423 

0.3487 

0.2 

1.770 

1  .084 

0.8846 

0  .  6975 

0 . 3 

2 . 655 

1.626 

1 . 327 

1.046 

0 . 4 

3.540 

2.167 

1.769 

1.395 

0 . 5 

4 .425 

2.709 

2 .212 

1  .744 

0  .  6 

5.310 

3.251 

2 . 654 

2.092 

0 . 7 

6.195 

3.793 

3 .096 

2.441 

00 

o 

7  .080 

4 .335 

3.538 

2.790 

0  .  9 

7 .965 

4  .877 

3 . 981 

3 .139 

1  . 

8  .  850 

5.419 

4 .423 

3.487 

1  . 1 

9 .735 

5 .961 

4.865 

3 .836 

1 . 2 

10  .  62 

6.502 

5.308 

4.185 

1 . 3 

11.50 

7  .044 

5 .750 

4 .534 

1  .  4 

12.39 

7.586 

6.192 

4.882 

Table  14.  Compton  Electron  FWHM  vs.  Incident  Gamma  Ray  FWHM  (y=500  keV) 


EEA 

FEP 

FEP 

FEP 

FEP 

FWHM 

FWHM 

FWHM 

FWHM 

FWHM 

(keV) 

ii 

CD 

(0=  .53) 

II 

CD 

(©=.  7) 

0 . 1 

0 .7233 

0.4490 

0 .3695 

0.2949 

0.2 

1  .447 

0.8981 

0.7391 

0.5898 

0 . 3 

2.170 

1 .347 

1.109 

0 .8846 

0 . 4 

2.893 

1.796 

1.478 

1.180 

0 . 5 

3.616 

2 .245 

1.848 

1 .474 

0  .  6 

4.340 

2 .694 

2 .217 

1.769 

0 . 7 

5 .063 

3 .143 

2 .587 

2.064 

00 

o 

5.786 

3.592 

2 . 956 

2 .359 

0  .  9 

6.510 

4.041 

3.326 

2 . 654 

1  . 

7 .233 

4.490 

3 .695 

2.949 

1  . 1 

7.956 

4.939 

4.065 

3 .244 

1 . 2 

8.679 

5.389 

4 .434 

3 .539 

1 . 3 

9.403 

5 .838 

4.804 

3 .833 

1  .  4 

10.13 

6.287 

5 .173 

4.128 
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Table  15.  Compton  Electron  FWHM  vs.  Incident  Gamma  Ray  FWHM  (y=600  keV) 


EEA 

FEP 

FEP 

FEP 

FEP 

FWHM 

FWHM 

FWHM 

FWHM 

FWHM 

(keV) 

ii 

CD 

(0  = ,  5  3) 

II 

CD 

(0=.7) 

0 . 1 

0  .  6156 

0 .3873 

0.3211 

0.2591 

0.2 

1 .231 

0 .7745 

0 . 6423 

0.5182 

0 . 3 

1  .847 

1.162 

0 .9634 

0 .7774 

0 . 4 

2.462 

1.549 

1 .285 

1 .036 

0 . 5 

3 .078 

1  .  936 

1.606 

1.296 

0  .  6 

3 .693 

2 . 324 

1  .  927 

1 .555 

0 . 7 

4.309 

2.711 

2 .248 

1.814 

00 

o 

4  .  924 

3 .098 

2.569 

2 . 073 

0  .  9 

5.540 

3.485 

2.890 

2 .332 

1  . 

6.156 

3 . 873 

3.211 

2.591 

1  . 1 

6.771 

4.260 

3.532 

2 .850 

1 . 2 

7 .387 

4  .  647 

3 .854 

3 .109 

1 . 3 

8 .002 

5 .034 

4 .175 

3 .369 

1  .  4 

8  .  618 

5 .422 

4.496 

3 . 628 

Table  16.  Compton  Electron  FWHM  vs.  Incident  Gamma  Ray  FWHM  (y=700  keV) 


EEA 

FEP 

FEP 

FEP 

FEP 

FWHM 

FWHM 

FWHM 

FWHM 

FWHM 

(keV) 

ii 

CD 

(0=  .53) 

II 

CD 

(©=.  7) 

0 . 1 

0 .5387 

0 .3432 

0.2866 

0 .2337 

0.2 

1  .077 

0.6864 

0 .5733 

0 .4674 

0 . 3 

1.616 

1 .030 

0 .8599 

0 .7010 

0 . 4 

2 .155 

1 .373 

1  .147 

0 . 9347 

0 . 5 

2 .693 

1.716 

1 .433 

1.168 

0  .  6 

3.232 

2 .059 

1.720 

1.402 

0 . 7 

3.771 

2.402 

2 .007 

1  .  636 

00 

o 

4.309 

2.746 

2.293 

1.869 

0  .  9 

4.848 

3.089 

2.580 

2.103 

1  . 

5 .387 

3 .432 

2.866 

2 .337 

1  . 1 

5 . 925 

3 .775 

3 .153 

2.571 

1 . 2 

6.464 

4.118 

3.440 

2.804 

1 . 3 

7 .003 

4.462 

3.726 

3 . 038 

1  .  4 

7.541 

4.805 

4 .013 

3 .272 
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Appendix  E:  CdTe  Detector  Systems  and  Supporting  Electronics 


A  total  of  seven  XR-lOOT-CdTe  detectors  with  associated  PX2T  power 
supply/shaping  amplifier  were  purchased  from  Amptek,  Inc.  System  7  was  the  original 
detector  system  used  by  Capt.  Williams.  This  system’s  PX2T  was  modified  for  a 
variable  high  voltage  bias  listed  below. 


Table  17.  List  of  CdTe  components 


CdTe  Detector  Systems  (EE A) 

System  # 

XR-100T 
Serial  Number 

PX2T 

Serial  Number 

Manufacture’s 
FWHM 
(122  keV) 

1 

03772 

2316 

1.145 

2 

03754 

2346 

1.319 

3 

03764 

2343 

1.24 

4 

03760 

2280 

1.427 

5 

03710 

2317 

1.239 

6 

03750 

2306 

1.084 

7 

03712 

2283 

1.181 

Additional  Electronics 

3125 

Canberra  Dual  HYPS 

1Z029517 
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Figure  65.  Photograph  of  XR-lOOT-CdTe  detector  with  housing  removed 


Figure  66.  Photograph  of  PX2T-CdTe  power  supply  and  shaping  amplifier 


99 


Appendix  F:  Attenuation  Calculations 


The  following  attenuation  calculations  were  conducted  to  analyze  the  effects  of 


the  incident  gamma  rays  entering  the  CdTe  through  the  aluminum  detector  housing 


versus  the  beryllium  window.  The  mass  attenuation  coefficients  and  densities  used 


below  were  obtained  from  XMuDat  (Nowotny,  1998). 

(* 

Beryllium  Data 
*) 

Gram 

Pee  =  1.848  - , 

(Centi  Meter)  3 

(Centi  Meter)  2 

P  MassBe  =  .06945  / 

Gram 

PBe  =  M  MassBe  PBe  ' 

BeAtten  =  (1  -  Convert  [10‘  Mil '  ] )  *  100  .  "%" 

0.325462  % 

(* 

Alluminum  Data 
*) 

Gram 

Pa l  =  2 . 699  - ; 

(Centi  Meter)  3 

(Centi  Meter)  2 

P  MassAl  =  .07566  -  ", 

Gram 

PA1  =  PMassAl  PA1  r 

AlAtten  =  (1  -  [1°-  1111  'Centi  *100.  "%" 


0.517341  % 
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Appendix  G:  Nal(Tl)  Detector  and  Supporting  Electronics 


Three  Nal(Tl)  detectors  with  integrated  photomultiplier  tubes  and  detachable 
preamplifiers  were  purchased  for  this  project  from  Saint-Gobain,  Inc.  To  operate  the 
PEA  the  high  voltage  supply  and  amplifier  listed  below  were  used. 


Table  18.  List  of  Nal(Tl)  components 


IN 

al(Tl)  Detectors  (PEA) 

Detector  # 

Detector 
Serial  Number 

Photomultiplier 

Tube 

Serial  Number 

Preamplifier 
PA- 14 

Serial  Number 

Manufacturer's 
Resolution 
(662  keV) 

1 

60004-00024-1 

75-4353 

70004-0063 

7.0% 

2 

60004-00025-1 

75-4317 

70004-0064 

7.2% 

3 

60004-00026-1 

75-4314 

70004-0065 

6.9% 

Additional  Electronics 

Model  9645 

Canberra  High  Voltage  Power  Supply 

08027948 

Model  672 

ORTEC  Spectroscopy  Amplifier 

329 

■ 


Figure  67.  Photograph  of  Nal(Tl)  detector 
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Figure  68.  Nal(Tl)  in  detector  housing  inside  lead  cave 


Appendix  H:  Coincidence  Electronics 

Two  methods  were  tested  for  determining  coincidence,  the  start-stop  method  and 
the  pulse  overlap  method.  Photographs  and  lists  of  electronics  follow. 


Figure  69.  Photograph  of  start-stop  coincidence  electronics 
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Table  19.  Electronics  used  for  start-stop  coincidence  technique 


TAC  Coincidence  Electronics 

Model 

Component 

Quantity 

Serial  Numbers 

Model  566 

ORTEC  Time-to-Amplitude  Converter/SCA 

1 

1547 

Model  551 

ORTEC  Timing  SCA 

1 

3880 

Model  579 

ORTEC  Fast-Filter  Amplifier 

2 

582/599 

Model  552 

ORTEC  Pulse-Shape  Analyzer/Timing  SCA 

2 

1166/1163 

Model  427A 

ORTEC  Delay  Amplifier 

3 

4369/3703/3441 

Figure  70.  Photograph  of  pulse  overlap  coincidence  electronics 
Table  20.  Electronics  used  for  pulse  overlap  technique 


LGSC  Coincidence  Electronics 

Model 

Component 

Quantity 

Serial  Numbers 

Model  409 

ORTEC  Linear  Gate  Slow  Coincidence 

1 

467 

Model  551 

ORTEC  Timing  SCA 

2 

3645/3957 

Model  416A 

ORTEC  Gate  and  Delay  Generator 

2 

004743/004397 

Model  427A 

ORTEC  Delay  Amplifier 

3 

4369/3703/3441 
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Appendix  I:  Initial  Collimator 


Two  initial  collimators  were  used  to  collimate  the  sample  source  to  the  EEA;  a 
manufactured  hexagonal  cross  section  collimator  and  a  locally  fabricated  circular  cross 
section  collimator.  The  first  collimator  was  used  during  Capt.  Williams’  research.  Both 
were  modified  to  include  a  source  holder  for  source  placement  reproducibility. 


Figure  71.  Photograph  of  hexagonal  cross  section  collimator 


Figure  72.  Photograph  of  circular  cross  section  collimator 
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Appendix  J:  Conical  Collimator 


The  conical  collimator  used  was  fabricated  by  Capt.  Williams  (Williams,  2003: 
122-123). 


Figure  73.  Photograph  of  conical  collimator  with  alignment  grid  in  center 
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Figure  74.  Photograph  of  conical  collimator  without  center  cone 
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Figure  75.  Photograph  of  conical  collimator  inside  cone  with  alignment  grid 


I 


Figure  76.  Dimensions  of  conical  collimator 
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Appendix  K:  Radioactive  Sources 


The  list  below  contains  all  sources  used  during  this  thesis. 


Table  21.  Radioactive  sources 


Source 

ID# 

Manf. 

Ref  Date 

Initial 

Activity 

Uncert. 

Active 

Diam. 

Cs-137 

T-083 

IPL* 

15-Jul-1998 

326.7  kBq 

±3.0  % 

3  mm 

Cs-137 

T-084 

IPL* 

15-Jul-1998 

375.2  kBq 

±3.0  % 

3  mm 

Co-57 

T-085 

IPL* 

15- Aug- 
1998 

3479  kBq 

±3.0  % 

3  mm 

Cs-137 

T-089 

IPL* 

1 -Aug- 1998 

370.7  kBq 

±3.1  % 

3  mm 

Multinuclide 

T-103 

IPL* 

15 -Mar- 
2001 

58.20  kBq 

±3.0  % 

5  mm 

Multinuclide 

T-105 

IPL* 

1 -Nov-2002 

190.3  kBq 

±3.1  % 

5  mm 

Na-22 

T-106 

IPL* 

15-Feb-2003 

32.87  kBq 

±3.1  % 

3  mm 

Na-22 

T-107 

IPL* 

15-Leb-2003 

32.55  kBq 

±3.1  % 

3  mm 

Multinuclide 

T-108 

IPL* 

15-Jul-2004 

37.67  kBq 

±3.4% 

5  mm 

Multinuclide 

T-109 

IPL* 

15-Jul-2004 

37.37  kBq 

±3.1  % 

5  mm 

Eu-152 

T-110 

IPL* 

15-Dec- 

2004 

375.2  kBq 

±3.0  % 

5  mm 

Na-22 

T-lll 

IPL* 

1 -Jan-2005 

1997  kBq 

±3.0  % 

3  mm 

*IPL-Isotope  Product  Laboratories 
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Appendix  L:  Sample  Calculation  of  Activity  Adjustment  for  Current  Date 


Standard  radioactive  decay  formula: 


Log[2.]  t 

A  [  t_  ,  Nuclide_]  =  Ao, Nuclide  ®  TNuclide 
y[t_,  Nuclide_]  =  BRNuclide  A  [ t ,  Nuclide] 


'  1  *10‘6  3.7  *1010  ) 

v  Micro  Curie  Second  ) 


Certificate  information  for  multinuclide  sample  T-105: 


(*  Multinuclide  T105  *) 
RefDateTios  =  {2002,  11,  1}; 

Ao , Am241  =  0 . 1522  Micro  Curie; 

TAm24i  =  Convert  [432 . 17  Year ,  Day]  ; 
BRAm241  =  -36; 

Ao , cdl 0 9  =  1.4  Micro  Curie; 

Tcdl09  =  4  62.6  Day; 

BRcdl09  =  .0363  ; 

Ao  ,co57  =  0.05457  Micro  Curie; 

TCo57  =  271.79  Day; 

BRco57  =  .  856  ; 

Ao,Tel23m  =  0 . 068  95  Micro  Curie; 
TTel23m  =  119.7  Day; 

BRTel23m  =  -84; 

Ao,cr51  =  1 . 836  Micro  Curie; 

TcrSI  =  27 . 706  Day; 

BRcrSI  =  .  0986  ; 

Ao , snll3  =  0 . 2530  Micro  Curie; 

Tsnll3  =  115.09  Day; 

BRsnll3  =  .648  9; 

Ao,sr85  =  0 . 3509  Micro  Curie; 

Tsr85  =  64.849  Day; 

BRsr85  =  .  984  ; 

Ao ,csl37  =  0 . 2365  Micro  Curie  ; 

Tcsl37  =  Convert  [30 . 17  Year ,  Day]  ; 
BRcs137  =  ■  851  ; 

Ao,y88  =  0 . 5066  Micro  Curie; 

Ty88  =  10  6.63  Day; 

BRy88  =  .  94  ; 

Ao ,co60  =  0 . 2840  Micro  Curie; 

Tco60  =  Convert  [5 . 272  Year,  Day]  ; 
BRco60  =  .  9986  ; 
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SampleDate  =  {2004,  7,  28}; 

DecayTime  = Day sBe tween [RefDateiios ,  SampleDate]  Day; 
A[DecayTime,  Am241] 

A  [DecayTime,  Cdl09] 

Y  [DecayTime,  Am241] 

Y  [DecayTime,  Cdl09] 

Am241  current  activity=  0.151776  Curie  Micro 
Cdl09  current  activity=  0.540644  Curie  Micro 

2021.66 

Am241  intensity=  Second 

726.139 

Cdl09  intensity=  Second 

Log  [2-]  t 

A  [  t_ ,  SId_,  Nuclide_]  =  A0,sid  ®  TNuclide  ; 

1  *  10~6  3.7  *  1010 

Y  [t_,  SId_,  Nuclide_]  =  BRh,v.h^  A[t,  Sid,  Nuclide]  - 

Micro  Curie  Second 

RefDateTo83  =  {1998,  07,  15}; 

RefDateT084  =  {1998,  07,  15}; 

RefDateTo89  =  {1998,  08,  01}; 

ao,t083  =  Convert  [326 . 7  Kilo  Becquerel ,  Micro  Curie]  ; 

Ao,t084  =  Convert  [375 . 2  Kilo  Becquerel ,  Micro  Curie]  ; 

Ao,t089  =  Convert  [370 . 7  Kilo  Becquerel ,  Micro  Curie]  ; 

Tcsi37  =  Convert  [30 . 17  Year ,  Day]; 

BRcsi37  =  .851; 

SampleDate =  {2004,  10,  28}; 

DecayTime [SId_]  = DaysBe tween [Ref Da teSid  ,  SampleDate]  Day; 

A [DecayTime [T083] ,  T083,  Csl37]  + 

A [DecayTime [T084] ,  T084,  Csl37]  + 

A [DecayTime [T089] ,  T089,  Csl37] 
r [DecayTime [T083] ,  T083,  Csl37]  + 

Y [DecayTime [T084] ,  T084,  Csl37]  + 

Y [DecayTime [T08 9] ,  T089,  Csl37] 

25.096  Curie  Micro 

790199. 

Second 
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Appendix  M:  Simulated  System  Computer  Code 


The  Mathematica  program  attached  below  was  written  to  simulate  the  Compton 
spectrometer  and  create  a  Compton  electron  energy  spectrum  as  seen  in  the  EEA.  The 
program  analyzes  the  effects  of  varying  geometric  parameters  on  the  resolution  and 
efficiency  of  the  Compton  spectrometer  system.  An  iterative  method  is  used  to  step 
through  a  finite  number  a  three-point  combinations,  see  Figure  77,  representing  a  gamma 
ray  emitted  from  an  incremental  source  area  then  Compton  scattered  in  an  incremental 
EEA  detector  volume  and  passing  through  an  incremental  area  at  the  bottom  of  the 
conical  collimator  where  it  is  absorbed  in  the  PEA.  If  the  photon  is  not  geometrically 
blocked  by  any  collimator  material  the  energy  of  the  resulting  Compton  electron  and  the 
probability  of  that  three-point  combination  occurring  are  determined. 


ill 


Check  points  along 
cone 

s Bottom  (Point  3) 


Figure  77.  Code  Simulation  geometry 
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Figure  78.  Example  of 
incremental  source  areas. 
Radius  pixel  set  to  3  and 
angle  pixel  set  to  8. 


Figure  79.  Example  of 
incremental  EEA 
volumes.  5x5x1  mm3  with 
proportional  pixel  set  to  1. 


Figure  80.  Example  of 
incremental  areas  of 
bottom  conical  collimator 
gap.  Radius  pixel  set  to  3 
and  angle  pixel  set  to  12. 


Z 


Figure  81.  Example  geometry  used  to  calculate  Compton  electron  energies  and 

Compton  scatter  probabilities. 


113 


The  energy  of  the  Compton  electron  is  calculated  using  the  Compton  relation 
described  in  Equation  2  and  an  estimated  scatter  angle  created  by  the  three-point 
combination.  Four  scatter  angles  are  calculated  for  each  three-point  combination  using 
the  center  point  of  the  incremental  source  area,  the  center  point  of  the  incremental 
detector  volume  and  the  four  comers  of  the  incremental  area  at  the  bottom  of  the  conical 
collimator  as  shown  in  Figure  8 1 .  The  average  of  these  four  scatter  angles  is  then  used  to 
calculate  the  Compton  electron  energy.  Similar  to  an  analog  to  digital  converter  (ADC), 
the  probability  of  this  combination  is  collected  in  an  “energy  bin”  based  on  the  Compton 
electron  energy  calculated  and  the  energy  versus  probability  spectrum  shown  in  Figure  82 
is  produced. 
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Figure  82.  Example  output  of  simulation  showing  Compton  electron  energy  versus 

probability  spectrum. 

The  probability  of  the  three-point  combination  occurring  is  the  product  of  the 
probability  that  the  incremental  area  of  the  isotropic  source  emits  a  gamma  ray  into  the 
top  area  of  the  incremental  detector  volume  and  the  probability  that  the  gamma  ray  will 
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Compton  scatter  from  the  incremental  detector  volume  into  the  incremental  area  at  the 
bottom  of  the  conical  collimator.  Treating  the  incremental  source  area  as  a  point, 
Equation  14  is  used  to  determine  the  solid  angle  subtended  by  the  incremental  detector 
volume.  To  calculate  the  scatter  probability  Equation  6  is  used  with  xl  and  x2 
determined  by  the  incremental  detector  volume  depth  and  p  calculated  using  the  Klein- 
Nishina  formula  shown  in  Equation  4.  The  formula  is  modified  for  the  sold  angle 
integration  over  d9  and  dtp  and  the  limits  of  the  integration  change  for  each  three-point 
combination.  The  scatter  angles  calculated  to  determine  the  Compton  electron  energy  are 
used  for  the  d9  limits  and  the  arc  size  of  the  conical  collimator  pixel  is  used  for  the  dtp 
limits. 


A  source  activity  and  collection  time  can  be  included  to  produce  the  typical 
Compton  electron  energy  versus  counts  spectrum  as  shown  in  Figure  83.  The 
incremental  source  areas  are  created  equal  therefore  the  activity  is  distributed  evenly  over 
the  entire  source  area. 
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Figure  83.  Example  output  of  simulation  showing  Compton  electron  energy  versus 

counts  spectrum. 
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«  Graphics’ ArgColors'; 

«  Graphics'Graphics3D'; 

«  Graphics'Shapes'; 

« UnearAlgebra’MatrixManipulation’; 

«  Miscellaneous’ Units’ 

Needs["Miscellaneous'RealOnly'"]; 

Off[General::"spelll"]; 

Off  [General::"spell"] ; 

ComptSpect  = 

Compile|{collz,  CollGap,  ZSource,  ZBotlnitColl,  {SourceAnglePix,  _Integpr},  {SourceRadPix,  -Integer}, 
{CollAnglePix,  -Integer},  {CollRadPix,  _Inte^r|, 

{PropPix,  -Integer},  {Checkpoints,  -Integer},  y,  DetWidthX,  DetWidthY,  DetWidthZ,  SourceRad, 
InitCollRad}, 

(*  ************  Initial  set  up  of  variables  ************  *) 

Act  =  (800*  1000)  *(3600*  l.)/((SourceAnglePix( SourceRadPix-  1))  +  1.); 

DetPixX  =  PropPixDetWidthX; 

DetPixY  =  PropPixDetWidthY; 

DetPixZ  =  PropPixDetWidthZ; 

DetWidthZ 

dDT  = - ; 

2.  DetPixZ 

ZTopInitColl  =  ZSource; 

Energy  Lower  =  60.; 

EnergyUpper  =  140.; 

Energy  BinSize  =  .1; 

EnergyBins  =  IntegerPart[(EnergyUpper  -  EnergyLower)/ Energy  BinSize]; 

Do[En[i]  =  0,  {i,  1,  EnergyBins}]; 

Do[NoProbEh[i]  =  0,  {i,  1,  EnergyBins}]; 

TopOuterConeCollRad  =  6.86; 

BotOuterConeCollRad  =  25.04; 

TopInnerConeCollRad  =  TopOuterConeCollRad  -  CollGap; 

BotlnnerConeCollRad  =  BotOuterConeCollRad  -  CollGap; 

ZTopConeColl  =  - 12.9  +  collz; 

ZBotConeColl  =  ZTopConeColl  -  31.0; 
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(*  ************  CfCtltC  SOUTCG  Radii  ************  *) 

Clear  [SourceR]; 

ClearfTenpRl; 

TempRfO]  =  0.; 

TempRfSourceRadPix]  =  SourceRad; 

TetnpRad  =  TablefO,  {i,  1,  2],  { j,  1,  2}]; 

Area  =  SourceRad2  /  ((SourceAnglePixfSourceRadPix-  1))  +  1.); 

TempRad=  Solve[Area=  TempRfi]2,  TempR[l]j; 

TeinpR[l]  =TempR[ll/.  TempRad[[2,  1]]; 

SourceRfll  =  0.; 

Do^TenpRad  =  Solve[Area  ==  (TempRfi]2  -  TempRfi -  l]2) / SourceAnglePix,  TempR[i]]; 
TempR[il  =  TeirpR[i]  /.  TempRad[[2,  1]]; 

SourceRfi]  =  Sqrt[|  —  (TempR[i]2  +  TempR[i  -  l]2)jj, 

{i,  2,  SourceRadPix-  l}j; 

SourceRfSourceRadPix]  =  Sqrt[|  —  (TenpR[SourceRadPix|2  +  TempRfSourceRadPix-  l]2)jj; 


(*  ************  Create  Source  Radii  ************  *) 

(*  ************  Create  Bottom  Collimator  Radii  ************  *) 

ClearfTenpRl; 

ClearfBottomCollRl ; 

TempRfO]  =  BotlnnerConeCollRad; 

TenpRfCollRadPix]  =  BotOuterConeCollRad; 

Area  = - (BotOuterConeCollRad2  -  BotlnnerConeCollRad2); 

CollRadPix 1  ’ 

Do[TenpRad  =  Solve[Area  ==  TempRfil2  -  TempRfi -  l]2,  TempRfi]]; 

TempRfil  =TenpR[i]/.  TenpRad[[2,  1]]; 

BottomCollRfi]  =  TempRfi], 

fi,  1,  CollRadPix-  1)]; 

BottomCollRfOl  =  BotlnnerConeCollRad; 

BottomCollRfCollRadPix]  =  BotOuterConeCollRad; 

(*  ************  Create  Bottom  Collimator  Radii  ************  *) 


(*  ************  Klein— Nishina  Formula  ************  *) 

Acd  =  112.411; 

Ate  =  127.60; 

ZcdTe  =  100;  (*  Electros  per  CdTe  *) 
pcdTe  =  6.06;  (*Density  from  Knoll*) 

6.022 14  xlO23 

NcdTe  = 


a[hv_]  = 


(Aqj  +  ATe) 

hv 

511’ 


r\— 15  . 


r0=  100  *  2.817940285*  10" 

(*  Note:  Z  is  the  number  of  electrons  *) 
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,  +  Cos[0]2  ) 

KN[0_,  Z_,  hv_]  =  Z  ro2  ( 1  /  ( 1  +  a[hv]  ( 1  -  Cos]#])))"  | - | 

l  2  I 

( 1  +  (ff[hv]2  ( 1  -  Cos[0])2)  /(( 1  +  Cos[0]2)  ( 1  +  a[hv]  ( 1  -  Cos]®))); 
(*  ************  Kleill— NistUlltl  FOlTOLlla  ************  *) 

(*  ************  END  Initial  set  up  of  variables  ************  *) 


(*  ************  Start  Calculations  ************  *) 

a=  0; 

(*  Detector  Z  *)Do[ 

(*  Detector  Y  *)Do[ 

(*  Detector  X  *)Do| 

(*  ************  Create  ConeCollTestData  Array  ************  *) 

(*  Collimator  Angle  *)Do| 

(*  Collimator  Radius  *)Do|^ 

(  r  2n  n'l 

XSlope  =  I  BottomCollR[CollRadIiidex]  Cos  - (CollAnglelndex  -  1)1  —  dpx; 

V  L  CollAnglePix  J) 

(  r  2n  A 

YSlope  =  I  BottomCollR[Col]RadIndex]  Sin - (CollAnglelndex-  1)  1 1  -  dpy; 

l  L  CollAnglePix  J) 

ZSlope  =  ZBotConeColl  -  dpz; 

Do[(*Test  points  in  channel*) 

ZCheck  =  ZTopConeColl  +  (ZBotConeColl  -  ZTopConeColl)  /  (QieckPoints  +  1)  k; 

RCheckMin  =  TopInnerConeCoURad  +  (BotlnnerConeCollRad  -  TopInnerConeCollRad)  /  (Checkpoints  +  1 )  k; 

RCheckMax=  TopOuterConeCollRad  +  (BotOuterConeCollRad-  TopOuterConeCollRad)  /  (Checkpoints  +  1)  k; 

ZCheck  -  dpz 

tt= - ; 

ZSlope 

(*Print[tt,"  ".RCheckMax2,"  ".RTest,"  ",  RCheckMin2];*) 

RTest  =  (dpx+  XSlopett)2  +  (dpy  +  YSlopett)2; 

If[ 

RTest  <  RCheckMax2  &&  RTest  >  RCheckMin2, 

(*  True  *) 

ConeTest  =  1, 

(*  False*) 

ConeTest  =  0., 

Print]] "Cone  Collimator  Test  Failed",  RCheckMax,  RCheckMin,  RTest|]; 

ConeTest  =  0. 

]; 

If  [ConeTest  ==  0,  Break]]], 

(k,  0,  CheckPointslj; 
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ConeCollTestData[ CollRadlndex  CollAnglelndex]  =  ConeTest; 

,  {CollRadlndex  0,  CollRadPix|]; 

,  {CollAnglelndex  1,  CollAnglePix}]; 

Do[ConeCollTestData[ CollRadlndex,  CollAnglePix+  11  =  ConeCollTestData[ CollRadlndex,  1], 

{CollRadlndex,  0,  CollRadPix)]; 

{ '■<  Create  ConeCollTestData  Array  ************  *) 

(*  *****  Test  ConeCollTestData  Array  Create  ConeCollTest  Array  *****  *) 

Do[If[ConeCollTestData[Col]RadIndex,  CollAnglelndex]  ==  1  && 

ConeCollTestData[CollRadIndex+  1,  CollAnglelndex]  ==  1  && 

ConeCollTestDatafCollRadlndex,  CollAnglelndex +  11  ==  1  && 

ConeCollTestData[CollRadIndex+  1,  CollAnglelndex  +  11  ==  1,  CollConeTest[ CollRadlndex,  CollAnglelndex]  =  1, 
CollConeTest[  CollRadlndex,  CollAnglelndex]  =  0]; 

( *Print  [  ConeT  est  [  CoURadIndex,CollAngleIndex|  ]  *) 

,  {CollRadlndex  0,  CollRadPix-  1],  {CollAnglelndex  1,  CollAnglePix)]; 

(*  *****  Test  ConeCollTestData  Array  Create  ConeCollTest  Array  *****  *) 

(*  Collintator*)Do| 

(*  Source  Radius  *)ixj 

If[CollConeTest[CollRadIndex  CollAnglelndex]  =0,  Break]]]; 

(*  Source  Angle  *) 

Do[ 

If[CollConeTest[CollRadIndex  CollAnglelndex]  ==0,  Break[]]; 

(*Initial  Collimator  Test  *) 

spx=  SourceR[SourceRadIndex]  Cos[SourceAng|; 

spy  =  SourceRfSourceRadlndex]  Sin[SourceAng|; 

XSlope  =  spx-  dpx 

YSlope  =  spy  -  dpy; 

ZSlope  =  ZSource  -  dpz; 

ZBotlnitColl-  dpz 

tt= - ; 

ZSlope 

RTest  =  (dpx+  XSlopett)2  +(dpy  +  YSlopett)2; 

If  RTest  <  InitCollRad2  , 


a  =  a+  1; 

(*  ************  Distance  traveled  in  CdTe  ************  *) 


ttDetTop  = 


0.  -  dpz 
ZSlope 


DT  =  Sqrt[(XSlopettDetTop)2  +  (YSlope  ttDetTop)2  +  dpz2 

(*  ************  Distance  traveled  in  CdTe  ************  *) 
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(*  ******  Solid  Angle:  From  Source  Subtended  by  Detector  ******  *) 
SolidAngle  = 

DetWidthY 

!  —  snx! !  !  dnv  + 

2.  DetPixX 


ArcTan 


ill  DetWidtbX  \ 

((  j  dpx+ 


+  -  _  _ 1  -  spy'1'1 

//  / 


(ZSource  -  dpz)  f  dpx  + 


SpX)[[dpy  2.  DetPixY 
DetWidthX 


2.  DetPixX 


\  X2  ( (  DetWidthY  t  t2 

|-SPX1  +iidpy+T5etPi^]-SPyi  +(ZS0UrCe 


-dpz) 


2)2 


ArcTan[j  jj  j'  dpx- 


DetWidtbX 
2.  DetPixX 


!  -  spx  dpy  - 


DetWidthY 
2.  DetPixY 


■  -  spy'1'1 , 
)  )  / 


If  r  DetWidthX 

(ZSource  -  dpz)  I  dpx - 

v  \\\  2.  DetPixX 


t  X2  ( (  DetWidthY  \  \2 

rspxi  +udpy+  i5etPi^]-spyj +(ZSource- 


A2)2  Ml 


dpz)zj”  !!  - 

;  )Y 


I  in  DetWidthX  t  \  ii  DetWidthY  t 

Hnv  4-  -  —  cn  v  Hn  v  —  -  !  —  cr\\/ 


ArcTanJj  dpx+ 


2.  DetPixX 


-spxjjjdpy- 


2.  DetPixY 


j-SPyjj/ 


III  DetWidthX 

(ZSource  -  dpz) !  dpx+ - 

v  VU  F  2.  DetPixX 


\  t2  ii  DetWidthY  t  \2 
|-SPX1  +  ||dpy  ZDetPbY^ J  _SpyJ  +(ZSource 


rf  if  f  DetWidthX  \  \(f  DetWidthY  t  \\ 

ArCTall[ll  i  dpX"  T^etPh^ )  -  SPXJ  UdPy  -  l5etPi^ J  -  HI/ 


(  fff  DetWidthX 

!  (ZSource  -  dpz) !  dpx - 

v  ill  F  2. DetPixX 


\  t2  ff  DetWidthY  \  \2 

]-spxi  ni^-T^etp ivrspy) +(ZSource 


-dpZ)f!!]+ 

;  j) 


-dpz)2l2]]l; 

j  )Y 


(*  ******  Solid  Angie:  From  Source  Subtended  by  Detector  ******  *) 

(*  **********  Calculate  scatter  angles  ***********  *) 

SP  =  {spx,  spy,  ZSource}; 

DP  =  {dpx,  dpy,  dpz}; 

CPI  =  |BottomCollR[CollRadIndex]  Cos[  ^  p-  (CollAngielndex  -  l)]; 

BottomColiR[CollRadIndex|  Sin[  ^  p-  (CollAngielndex-  1)],  ZBotConeColl]; 


CP2  =  |BottomCollR[CollRadIndex  +  11  Cos 


2  n 


—  (CollAngielndex-  1)1, 

Jix 


-  CollAnglePix 

r  2  n  , 

BottomColiR[CollRadIndex+  11  Sin - (CollAngielndex-  1)1,  ZBotConeColl}; 

L  CollAnglePix  J  ’ 

CP3  =  |BottomCollR[CollRadIndex]  Cos[  ^  p-  (CollAngielndex)], 


BottomColiR[CollRadIndex]  Sin} 


In 


CollAnglePix 


(CollAngielndex)],  ZBotConeColl}; 
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CP4  =  |BottomCoUR[CollRadIndex+  1]  Cos[ 


2  n 


—  (CollAnglelndexjJ, 


-  CollAnglePix 

BottoiiiCollR[CoIIRadIndex+  1]  Sin[  ^  p-  (CollAngleIndex)J,  ZBotConeCollJ; 


ScatAnglel  =(n  -  ArcCos[ 

Dot[SP  -  DP,  CPI  -  DP]  /  (Sqrt[Dot[SP  -  DP,  SP  -  DP]]  Sqrt[Dot[CPl  -  DP,  CPI  -  DP]])]); 
ScatAngle2  =  (n  -  ArcCos[ 

DotfSP  -  DP,  CP2  -  DP]  /  (Sqrt[Dot[SP  -  DP,  SP  -  DP]]  Sqrt[Dot[CP2  -  DP,  CP2  -  DP]])]); 
ScatAng]e3  =  (n  -  ArcCos[ 

DotfSP  -  DP,  CP3  -  DP]  /  (SqrtfDotfSP  -  DP,  SP  -  DP]]  Sqrt[Dot[CP3  -  DP,  CP3  -  DP]])]); 
ScatAngle4  =  (n  -  ArcCosf 

DotfSP  -  DP,  CP4  -  DP]  /  (SqrtfDotfSP  -  DP,  SP  -  DP]]  Sqrt[Dot[CP4  -  DP,  CP4  -  DP]])]); 

(*  **********  Calculate  scatter  angles  ***********  *) 

(*  **********  Calculate  Compton  Cross  Sections  ***********  *) 


lft  =  Mntegute|^NQiTeKN[0,  ZojTe,  662.]  *  Sin[0],  (6,  ScatAngJe3,  ScatAngle4],  | </>,  0, 


2  71 


CollAngJePix^’ 


rt  =  Nlntegrate^NcdTeKNfO,  Zqit&  662.]  *  Sinfff],  [6,  ScatAnglel,  ScatAngle2],  j 0,  0, 


2  n 

CollAnglePix 


}]  = 


lft  +  rt 
2. 


(*  **********  Calculate  Compton  Cross  Sections  ***********  *) 

(*  **********  Calculate  Compton  Scatter  Probability  ***********  *) 


ScatProb  =  e1* PCdTe '  10.  ,-e- 
(*  **********  Calculate  Compton  Scatter  Probability  ***********  *) 

(*  **********  Calculate  Compton  Electron  Energy  ***********  *) 

1 

Scatter  Angle  =  —  (ScatAnglel  +  ScatAngJe2  +  ScatAng)e3  +  Scat  Angle4) ; 

/y 

( 1  +  - —  ( 1  -  Cos[ScatterAngle])J; 


I  DT-dDT\ 


DT+dDT  \ 


I  Ul+cl 

_-!•  PCdTe  (^a 


Enlndex  =  IntegsrPart^ 


511. 

Energy  Deposited  -  Energy  Lower 
Energy  BinSize 


EnfEnlndex]  =  EnfEnlndex]  +  (ScatProb  *  SolidAngle); 

NoProbEnf  Enlndex]  =  NoProbEiifEnlndex]  +  1.; 

( *  **********  Calculate  Compton  Electron  Energy  ***********  * 
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If[SourceRadIndex==  1,  Break[H,  jSourceAng 

— — — 11; 

SourceAnglePix  ’ 

,  {SourceRadlndex,  1,  SourceRadPix)  ; 


n  ( 

—  ,12.; 


SourceAnglePix  l 


,  {CollRadlndex,  0,  CoURadPix-  1},  {ColiAngJelndex,  1,  ColLAnglePix)|; 


l. 

DetWidthX  DetWidthX 

DetWidthX 

DetWidthX  ( 

|dpx, 

2.  '  2.DetPixX  ’ 

2. 

DetPixX 

{dpy>  - 

DetWidthY  DetWidthY 

DetWidthY 

DetWidthY ) , 

2.  '  2.  DetPixY  ’ 

2. 

DetPixY  A’ 

Priiit[dpz]; 


.  {dpz. 


DetWidthZ 
2.  DetPixZ 


-DetWidthZ, 


DetWidthZ 

DetPixZ 


n  \ 

- I, 

SourceAnglePix ) 
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Appendix  N:  CdTe  Spectra 


191 1 
1764 
1617 
1470 
1323 
1176 
1029 
882 
735 
588 
441 
294 
147 


. J 

J 

Figure  84.  57Co  spectrum  using  CdTe 
detector  (RTD  off,  Position  1,  Figure 
20) 


Figure  87.  57  Co  spectrum  using  CdTe 
detctor  (RTD  on,  Position  2) 
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Figure  85.  57Co  spectrum  using  CdTe 
detector  (RTD  off,  Position  2) 
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Figure  88.  57Co  spectrum  using  CdTe 
detector  provided  by  Amptek 
(Amptek,  2004:  25) 
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Figure  86.  57Co  spectrum  using  CdTe 

detector  (RTD  on,  Position  1)  Figure  89.  137Cs  spectrum  using  CdTe 

detector  (RTD  off,  Position  1) 
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Figure  90.  137Cs  spectrum  using  CdTe 
detector  (RTD  off,  Position  2) 
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Figure  93.  137Cs  spectrum  using  CdTe 
detector  provided  by  Amptek 
(Amptek,  2004:  28) 


Figure  91.  137Cs  spectrum  using  CdTe 
detector  (RTD  on,  Position  1) 
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Figure  92.  137Cs  spectrum  using  CdTe 
detector  (RTD  on,  Position  2) 


125 


Appendix  O:  Nal(Tl)  Spectra 
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Figure  94.  Multinuclide  (T108)  spectrum  using  Nal(Tl)  detector 


126 


927 


Figure  95.  Nal(Tl)  detector  characterization  spectrum;  57Co 


Figure  96.  Nal(Tl)  detector  characterization  spectrum;  137Cs 
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Appendix  P:  CdTe  Detector  Diagram 


Figure  97.  Compton  Detector  System  Geometry 
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Figure  98.  Photograph  of  Compton  spectrometer  in  situ 
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